LASER  BIOEFFECTS  RESULTING  FROM  NON-LINEAR 
INTERACTIONS  OF  ULTRASHORT  PULSES  WITH 
BIOLOGICAL  SYSTEMS 


Final  Report 

AFOSR  Grant  F49620-01-1-0211 


Randolph  D.  Glickman,  Ph.D.,  P.I. 

Department  of  Ophthalmology 
University  of  Texas  Health  Science  Center  at  San  Antonio 
San  Antonio,  Texas 


Grant  Period;  1  April  2001  to  14  July  2004 


"SSffissr  20050427  085 

Distribution  Unlimited 


i 


TABLE  OF  CONTENTS 


Summary  .  1 

Section  One 

“Multiphoton  Absorption  and  Ultrashort  Pulse  Laser  Bioeffects  in  the  Retinal  Pigment 
Epithelium” . 3 

Section  Two 

“Intracellular  Signaling  Mechanisms  Responsive  to  Laser-Induced  Photochemical  and  Thermal 
Stress” . 21 

Section  Three 

“Nonin  vasive  Thermography  of  Laser-Induced  Hyperthermia  Using  Magnetic  Resonance”  ...  38 
Appendix  of  Supplemental  Information  .  60 


li 


Laser  Bioeffects  Resulting  from  Non-Linear  Interactions 
of  Ultrashort  Pulses  with  Biological  Systems 

Final  Report 

AFOSR  Grant  F49620-01 -1-021 1 

Randolph  D.  Glickman,  Ph.D. 

Department  of  Ophthalmology 
University  of  Texas  Health  Science  Center  at  San  Antonio 
San  Antonio,  Texas 


Summary 


The  original  goal  of  this  project  was  to  determine  the  role  of  non-linear  interactions  underlying  the 
bioeffects  induced  by  ultrashort  pulse  laser  pulses.  As  initially  conceived,  this  line  of  investigation 
was  to  be  principally  directed  at  understanding  the  contribution  of  multiphoton  absorption.  This 
indeed  was  a  major  focus  of  the  research  project,  but  for  various  reasons  the  scope  of  the  work  was 
expanded  to  include  identification  of  the  intracellular  mechanisms  that  determine  the  cellular 
response  to  the  absorption  of  optical  radiation,  and  to  develop  and  implement  a  non-invasive  means 
for  measuring  the  thermal  gradients  induced  by  the  absorption  of  laser  radiation  in  tissue.  Both  of 
these  ancillary  projects  were  successful  in  that  (1)  the  transcription  factor  NF-B  was  found  to  be 
activated  by  visible  laser  exposure  in  a  way  that  appeared  to  be  dependent  on  the  absorption  of  laser 
energy  in  the  melanin  granules  of  the  retinal  pigment  epithelial  cell,  and  (2)  by  exploiting  the 
temperature-dependent  nature  of  the  proton  resonance  frequency  (PRF),  magnetic  resonance 
thermography  was  successfully  used  to  measure  temperature  gradients  induced  in  tissue  phantoms 
during  laser  exposure,  and  these  gradients  closely  followed  the  spatial  distributions  predicted  by 
classical  heat  diffusion  theory. 


In  the  area  of  non-linear  tissue  interaction  mechanisms,  laser-induced  breakdown  has  the 
lowest  energy  threshold  in  tSe  femtosecond  domain,  and  is  responsible  Tdrf)fo(iuctron'bf  threshold 
ocular  lesions.  It  has  been  proposed  that  multiphoton  absorption  may  also  contribute  to  ultrashort- 
pulse  tissue  damage,  based  on  the  observation  that  33  fs,  8 1 0  nm  pulse  laser  exposures  caused  more 
DNA  breakage  in  cultured,  primary  RPE  cells,  compared  to  C  W  laser  exposures  delivering  the  same 
average  power.  DNA  breakage  was  assessed  with  single  cell  gel  electrophoresis,  i.e.  the  “comet” 
assay.  Subsequent  studies,  demonstrating  two-photon  excitation  of  fluorescence  in  isolated  RPE 
melanosomes,  appeared  to  support  the  role  of  multiphoton  absorption,  but  mainly  at  suprathreshold 
irradiance.  Additional  experiments,  however,  did  not  show  a  consistent  difference  in  the  DNA 
strand  breakage  produced  by  ultrashort  and  CW  threshold  exposures.  DNA  damage  appeared  to  be 
dependent  on  the  amount  of  melanin  pigmentation  in  the  cells,  rather  than  the  pulsewidth  of  the 
laser;  current  studies  have  found  that,  at  threshold,  CW  and  ultrashort  pulse  laser  exposures 
produced  almost  identical  amounts  of  DNA  breakage.  A  theoretical  analysis  suggested  that  the 
number  of  photons  delivered  to  the  RPE  melanosome  during  a  single  250  fsec  pulse  at  the  ED}0 
irradiance  is  insufficient  to  produce  multiphoton  excitation.  This  result  appeared  to  exclude  the 
melanosome  as  a  locus  for  two-  or  three-photon  excitation;  however,  a  structure  with  a  larger 
effective  absorption  cross-section  than  the  melanosome  may  interact  with  the  laser  pulses.  One 
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possibility  is  that  the  nuclear  chromatin  acts  as  a  unit  absorber  of  photons  resulting  in  DNA  damage, 
but  the  near  equivalence  of  the  comet  assay  results,  in  RPE  cells  exposed  to  ultrashort  and  C  W  laser, 
indicates  that  multiphoton  absorption  is  not  a  major  contributor  to  the  ultrashort  pulse  laser  damage 
threshold  in  the  near  infrared  (but  is  likely  to  contribute  to  suprathreshold  tissue  damage). 

The  response  of  biological  cells  to  laser  radiation  is  the  result  of  a  complex  series  of 
interactions.  The  cells  of  the  retinal  pigment  epithelium  (RPE)  are  subject  to  photo-oxidative  stress 
arising  from  the  interaction  of  incident  light  with  iipofuscin,  melanin,  and  other  pigment  granules 
in  the  RPE  cytoplasm.  Specific  genotypic  responses  to  these  stressors  are  controlled  by  transcription 
factors,  such  as  NF-kB  (RelA/p50  dimer).  The  effect  of  CW  laser  exposures  on  NF-kB  nuclear 
translocation  have  been  studied  in  a  line  of  human-derived  RPE  cells  (hTERT-RPE)  that  develop 
melanin  pigmentation  in  culture.  The  cells  were  exposed  to  the  C  W  emission  of  an  Argon-ion  laser 
for  10  m  at  0.5  W/cm2,  a  range  previously  shown  to  produce  oxidation  of  cellular  proteins,  DNA, 
and  antioxidants.  NF-kB  dimer  was  measured  in  nuclear  extracts  by  an  electrophoretic  mobility 
shift  assay.  NF-kB  nuclear  translocation  exhibited  a  modest,  early  peak  at  1  h,  and  a  larger,  late 
peak  at  24  h.  NF-kB  activation  could  be  reduced  only  by  some  antioxidants;  for  example,  20  mM 
N-acetyl-L-cysteine  or  100  pM  pyrrolidine  dithiocarbamate  were  ineffective,  while  500  pM 
ascorbic  acid  was  highly  effective.  These  results  indicate  that  interaction  of  the  laser  with  the  RPE 
melanin  granules  is  a  likely  source  of  oxidative  reactions,  and  the  induction  of  photoxidative  stress 
activates  NF-kB,  but  it  remains  to  be  determined  if  NF-kB  is  pro-  or  anti-apoptotic  in  the  RPE  cell. 


An  important  component  of  laser-tissue  interaction  is  the  production  of  thermal  stress  in  the 
radiated  tissue.  The  possibility  to  induce  selective  hyperthermia  in  a  target  tissue  or  organ  is  of  great 
interest  for  the  treatment  of  cancer  and  other  diseases,  however,  the  therapy  is  currently  limited 
because  the  temperature  required  for  optimal  tissue  response  is  unknown,  and  the  ability  to  perform 
accurate  thermography  in  the  radiated  tissue  may  be  limited.  To  overcome  this  limitation,  we 
implemented  a  non-invasive  thermographic  method  and  validated  it  by  measuring  laser-induced 
temperature  gradients  in  an  ocular  tissue  phantom.  Magnetic  resonance  thermography  (MRT)  was 
used  as  a  non-invasive  method  to  determine  the  temperature  distribution  inside  the  phantom  during 
exposure  to  a  continuous  wave  diode  laser  at  806  nm  wavelength  with  1  watt  maximum  output.  The 
laser  beam  had  a  qunsi-gaussian  profile,  with  a  radius  of  0. 8-2.4  mm  at  the  target.  High  quality 
temperature  images  were  obtained  from  temperature-dependent  phase  shifts  in  the  proton  resonance 
frequency  with  a  resolution  of  ldeg  C  or  better,  using  a  2T  magnet.  A  phantom  with  a  layer  of 
bovine  RPE  melanin  of  1.5  mm  thickness  was  used  to  determine  the  spatial  resolution  of  the  MRT 
measurements.  Three  dimensional  temperature  maps  were  also  constructed  showing  a  spatial 
resolution  of  0.25  mm  in  all  directions.  The  heat  distribution  depended  on  the  laser  parameters,  as 
well  as  the  orientation  of  the  melanin  layer  with  respect  to  the  incident  laser  beam.  The  temperature 
profiles  determined  by  MRT  closely  followed  predictions  of  a  heat  diffusion  model,  based  on  the 
optical  properties  of  infrared  light  in  melanin.  These  results  support  the  use  of  MRT  to  map 
temperature  gradients  in  a  small  organ  such  as  the  eye,  as  well  as  to  validate  models  of  laser-induced 
thermal  tissue  damage. 
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SECTION  ONE 


Multiphoton  Absorption  and  Ultrashort  Pulse  Laser  Bioeffects  in  the  Retinal  Pigment 

Epithelium 


INTRODUCTION 

Multiphoton  absorption 

The  possibility  of  multiphoton  absorption  was  first  hypothesized  by  Goppert-Mayer  in  1931, 
who  proposed  a  mechanism  by  which  photons  absorbed  within  a  sufficiently  short  time  window, 
i.e.  <10'18  sec,  could  cumulatively  add  their  energy  to  a  chromophore,  exciting  it  to  an  energy 
level  approximately  twice  as  great  as  that  of  the  individual  photons.1  Due  to  the  two-photon 
absorption  mechanism,  relatively  low  energy  photons  delivered  in  a  short  pulse  could  produce  an 
effect  equivalent  to  the  absorption  of  a  single,  more  energetic  photon.  The  practical  significance 
of  the  two-photon  excitation  mechanism  was  not  realized  until  ultrashort  pulse  lasers  were 
available.2-4 

Relevance  of  multiphoton  absorption  for  damage  mechanisms  involved  in  ultrashort  pulse  laser 
tissue  interactions 

In  previous  work  on  the  interaction  of  pulsed,  infrared  laser  sources  on  retinal  pigment  epithelial 
(RPE)  cells,  we  demonstrated  that  effects  normally  requiring  visible  or  near  UV  (UVA)  light 
exposure,  such  as  excitation  of  melanin  fluorescence5  and  DNA  breakage,6,7  could  be  produced 
by  ultrashort  infrared  (810  nm)  laser  pulses.  It  was  proposed  that  these  effects  were  due  to 
muitiphoton  absorption  within  the  melanin  bf  the  RPE  ceils.  Other  authors  have  also  reported 
two-photon  excitation  of  melanin  fluorescence.8,9  Two-photon  excitation  effectively 
upconverted  the  infrared  photons  to  UVA- VIS  wavelengths,  which  then  activated  mechanisms 
that  induced  single  and  double  strand  DNA  breaks.7,10  Although  an  increased  number  of  DNA 
strand  breaks  were  demonstrated  at  threshold  radiant  exposures,  the  differences  observed 
between  the  effects  of  continuous  wave  exposures,  which  would  not  be  expected  to  induce 
multiphoton  effects,  and  ultrashort  pulse  exposures,  which  might  produce  bioeffects  due  to 
multiphoton  absorption,  have  been  difficult  to  document  consistently.  In  this  report,  we  present 
our  experimental  data  on  the  ability  of  ultrashort  laser  pulses  to  break  genomic  DNA,  as  well  as  a 
theoretical  prediction  of  the  number  of  photons  likely  to  be  absorbed  in  cellular  targets  during 
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short  laser  pulses.  This  analysis  indicates  that  multiphoton  absorption  in  the  RPE  melanosome  is 
not  likely  to  contribute  to  threshold  damage,  although  it  probably  contributes  to  suprathreshold 
cellular  damage  due  to  induction  of  photo-oxidative  stress,  for  example,  as  shown  by  Denton  et 
al." 

METHODS 

The  hTERT-RPE  cell  line  and  comet  assay  for  DNA  breakage 

Experiments  were  performed  on  a  line  of  human-derived  RPE  cells  provided  by  the  Geron  Corp. 
(Menlo  Park,  CA).  These  cells,  designated  hTERT-RPE,  have  been  transfected  with  the  gene 
for  human  telomerase  (hTERT)  and  do  not  become  senescent  in  culture;  however,  in  other 
respects  the  cells  exhibit  a  normal  phenotype.12,13  The  hTERT-RPE  cells  produce  melanin 
pigmentation  in  culture,14  and  most  of  the  cells  used  in  these  experiments  exhibited  a  degree  of 
pigmentation  proportional  to  their  time  in  culture.  The  hTERT-RPE  cells  were  grown  for  at  least 
2  weeks  in  24-well,  plastic  culture  plates,  in  standard  Dulbecco’s  Modified  Eagle’s  Medium 
(DMEM)  containing  10%  serum.  During  the  laser  exposures,  the  overlying  media  was  drained 
from  the  wells  and  replaced  with  Dulbecco’s  Phosphate  Buffered  Saline  (DBPS).  Further  details 
of  the  procedures  for  handling  the  hTERT-RPE  cells  were  described  elsewhere.7 

Assessment  of  DNA  breakage  by  comet  assay 

The  use  of  the  comet  assay  for  detecting  DNA  strand  breaks  in  RPE  cells  was  reported  by  Hall  et 
al.6  Following  electrophoresis  of  the  lysed  RPE  cells  in  agarose  gels,  DNA  fragments  were 
stained  with  Syber  Yeiiov?;  a  fiubrescent  nucleic  acid  stdin  (Moiecular  Probes,  Eugene,  OR): 
The  fluorescent  DNA  patterns  (comets)  were  imaged  with  a  fluorescence  microscope  equipped 
with  a  CCD  imaging  system  (Optronics  DEI-470),  and  the  images  processed  with  a  public 
domain,  comet  assay  analysis  macro  for  the  NIH/Scion-Image  software15,  in  order  to  extract 
comet  tail  length,  tail  moment  (product  of  length  and  intensity),  and  comet  area.  Statistical 
analysis  of  the  results  was  done  using  the  Bonferroni  multiple  comparison  test  in  the  ProStat 
statistical  package  (Poly  Software  International,  Pearl  River,  NY). 

A  modification  was  made  to  the  analysis  algorithm  to  address  the  problem  that  high-intensity 
fluorescent  figures,  such  as  those  from  control  samples,  were  scored  by  the  algorithm  as  having 
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very  large  comet  tail  lengths  and  moments.  The  algorithm  used  in  our  earlier  experiments7 
sometimes  reported  larger  comet  tails  and  moments  for  the  cells  in  the  control  conditions  than 
for  cells  in  the  experimental  conditions,  although  by  eye,  the  control  comets  appeared  to  be 
“rounder”  and  less  comet-like.  The  original  comet  analysis  algorithm15  operated  by  thresholding 
the  comet  image,  then  measuring  the  length  and  width  of  the  border  of  the  threshold  image.  We 
originally  used  the  same  absolute  threshold  for  all  comets  in  a  given  experiment,  to  avoid  using  a 
subjective  setting  for  the  analysis.  This  introduced  an  artifact,  however,  by  over-emphasizing  the 
size  of  bright  comets,  when  the  threshold  was  set  to  detect  the  outline  of  less  intense  comets. 
This  was  due  to  the  fact  that  the  comets  in  the  control  conditions  tended  to  be  very  intense  (the 
DNA  did  not  migrate  in  the  gel,  but  remained  near  the  location  of  the  cell  nucleus),  while  the 
comets  in  experimental  conditions  tended  to  be  relatively  less  intense  (the  DNA  with  an 
increased  number  of  breaks  migrated  through  the  gel  during  electrophoresis,  and  therefore  was 
more  dilute).  To  minimize  this  problem,  the  analysis  algorithm  was  modified  so  that  the 
threshold  for  a  given  comet  was  set  to  65%  of  the  maximum  white  level  of  that  particular  image. 
In  effect,  the  threshold  used  to  determine  the  comet  length  was  determined  by  a  fixed  ratio  based 
on  the  dynamic  range  of  the  comet  image,  not  at  an  absolute  level.  The  65%  level  was  derived 
empirically,  i.e.  it  was  the  threshold  setting  that  did  not  over-emphasize  the  size  of  the  control 
comets,  while  still  detecting  the  comet  tails  in  the  experimental  conditions.  This  modification 
was  justified  because  the  absolute  white  level  and  contrast  range  varied  among  different  cell 
samples. 

Melaiiosoms  preparation  and  hTERT-RPE  cell  mslanosoms  loading 

RPE  cells  were  isolated  from  bovine  eyes  purchased  from  AnTech  (Tyler,  TX).  Following 
sonication  to  release  the  cells'  contents,  the  melanosomes  were  isolated  by  successive 
centrifugation  using  methods  previously  reported.16  The  melanosomes  were  stored  frozen  in 
0.25  M  sucrose  at  -20°  C  until  use.  The  stock  melanosome  suspension  contained  approximately 
6  x  109  (±20%)  granules/ml,  as  determined  by  haemocytometer  counts.  Melanosomes  were 
layered  onto  confluent  cultures  of  hTERT-RPE  24  h  prior  to  laser  exposures,  following  the 
method  of  Denton  et  al.11  Excess  melanosomes  were  removed  prior  to  laser  exposure  by  rinsing 
the  cells  in  normal  DMEM. 


5 


Laser  sources  and  exposures 

The  laser  source  used  in  the  hTERT-RPE  cell  exposure  experiments  was  a  Coherent  (Santa 
Clara,  CA)  MIRA  900F  Ti:Sapphire  regenerative  amplifier  system,  modified  from  the  system 
originally  described  by  Noojin  et  al.17  This  system  could  be  operated  to  produce  either  CW 
output  at  810  nm,  or  mode-locked  pulses.  In  mode-locked  (ML)  operation,  it  was  anticipated 
that  the  laser  produced  an  80  MHz  output  train  of  ~48-fsec  pulses  with  a  40  nm  bandwidth 
centered  at  810  nm.  Subsequent  beam  analysis,  however,  revealed  that  the  pulse  widths 
produced  were  actually  between  200  and  300  fsec.  This  discrepancy  may  have  played  a  critical 
role  in  the  determination  of  the  threshold  damage  mechanisms,  as  will  be  discussed  later  in  this 
section. 

An  external  lens  was  used  to  focus  the  810  nm  beam  onto  the  sample  plane.  The  beam  diameter 
at  1/e  was  0.210  mm  at  the  RPE  cell  monolayer  on  the  culture  plate.  The  laser  was  delivered  to 
the  cells  in  an  exposure  of  0.25  s  duration,  and  the  average  power  of  both  the  CW  and  the  train  of 
mode-locked,  ultrashort  pulses  was  set  to  either  80  or  160  W/cm2  depending  on  the  experiment. 
These  values  were  estimated  as  the  irradiance  equivalent  to  the  corneal  irradiance  corresponding 
to  0.5  x  ED50  and  1  x  ED50,  respectively,  for  a  retinal  lesion  in  the  intact  (human)  eye.  The 
average  power  of  both  laser  sources  was  measured  with  a  Molectron  (Portland,  OR)  EPM  1000 
meter  equipped  with  a  PowerMax  PM- 10  thermopile  detector.  Just  prior  to  the  onset  of  the  laser 
exposures,  the  culture  well  was  drained  of  excess  culture  medium,  and  replaced  with  100  pi  of 
DPBS,  To  expose  as  many  cells  in  each  culture  well  as  possible,  th,e  laser  was  scanned  over. the 
well,  delivering  the  0.25  s  exposures  in  a  spiral  pattern.  Approximately  1,261  exposures  were 
required  to  expose  each  well. 

The  cells  were  exposed  to  the  810  nm  CW  and  ML  laser  as  follows:  (1)  controls  received  no 
laser  exposure,  but  were  handled  in  a  similar  fashion  to  the  laser-exposed  cells;  (2)  exposed  to 
the  810  nm  ML  pulses  at  the  estimated  ED50  exposure  required  to  produce  a  retinal  lesion  (~160 
W/cm2)  (note:  the  pulses  were  originally  thought  to  be  48-fsec  pulsewidth,  but  were  actually 
~250-fs  at  80  MHz  in  a  train  of  0.25  sec  duration),  (3)  exposed  to  the  810  nm  ML  pulses  at  0.5  x 
ED50  (~80  W/cm2),  (4)  exposed  to  810  nm  laser  in  CW  mode  at  the  ED50  (0.25  sec  at  160 
W/cm2),  or  (5)  exposed  to  810  nm  CW  mode  laser  at  0.5  x  ED50  (0.25  sec  at  80  W/cm2). 
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RESULTS 


Evidence  for  multiphoton  absorption  in  the  RPE  cells 

Evidence  for  multiphoton  absorption  produced  by  the  ultrashort  pulse  laser  exposures  was 
provided  by  the  observation  of  melanin  fluorescence  excited  by  810-nm  ML  pulses  in 
melanosomes  isolated  from  RPE  cells.5  The  fluorescence  observed  in  these  experiments 
exhibited  a  similar  wavelength  maximum  at  approximately  525  nm  as  that  elicited  by  single¬ 
photon  excitation  with  the  406  nm  output  of  a  Kr-ion  CW  laser.  The  power  dependence  of  the 
two-photon  excited  fluorescence,  however,  closely  followed  a  second  order  function,  as  opposed 
to  single-photon  excited  fluorescence,  which  was  linearly  dependent  on  the  excitation  power 
(Figure  1).  This  observation  was  consistent  with  two-photon  excitation  of  melanin  fluorescence 
by  the  ultrashort  pulse,  810  nm  laser. 


Figure  1:  Power  dependence  of  fluorescence  in  isolated  RPE  melanosomes 
DNA  damage  by  ultrashort  and  CW  infrared  laser  exposure 

These  experiments  compared  the  amount  of  DNA  damage  produced  in  cultured  hTERT-RPE 
cells  by  exposures  to  810  nm  CW  emissions  and  to  ultrashort  mode- locked  (ML)  laser  pulses 
delivering  equivalent  average  power  to  the  cells.  Two  exposure  conditions  were  examined:  one- 
half  the  ED5o  fluence  required  for  retinal  lesion  threshold  (80  W/cm2),  and  threshold  (ED5o) 
fluence  (160  W/cm2).  The  experiments  were  conducted  with  hERT-RPE  cells  not  supplemented 
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with  melanosomes.  These  cells  spontaneously  produce  melanin  pigment  in  their  cytoplasm,14 
but  were  only  lightly  pigmented  at  the  time  these  experiments  were  carried  out.  The  effects  of 
the  laser  exposures  on  the  cells  were  assessed  by  determining  single  and  double  DNA  strand 
breaks  released  by  alkaline  hydrolysis  of  the  cells,  using  the  comet  assay  procedure  described 
above  (“Assessment  of  DNA  breakage  by  comet  assay”).  The  results  of  the  0.5  x  EDS0 
exposures  are  shown  in  Figure  2  A,  and  the  results  of  the  ED50  exposures  are  shown  in  Figure  2B. 
At  0.5  x  ED50,  cells  exposed  to  ML  laser  pulses  had  slightly  longer  comet  tail  lengths  than  did 
the  CW-exposed  cells  (Figure  2A).  The  difference  was  significant  (p<.001,  Bonferroni  multiple 
comparison  test).  The  comet  tail  lengths  of  the  CW-exposed  cells  were  not  significantly 
different  from  those  of  the  control  cells.  At  the  ED50  (Figure  2B),  the  comet  tail  lengths  in  both 
the  CW  and  ML  groups  were  significantly  longer  than  in  the  control  cells  (p<005),  and  the  tail 
lengths  of  the  ML-exposed  cells  were  longer  than  those  of  the  CW-exposed  cells  (p<.02), 
suggesting  that  the  ML  pulses  were  more  effective  in  producing  DNA  strand  breaks.  The  data, 
however,  were  somewhat  inconsistent,  because  the  comet  tail  moments,  a  measure  of  the  density 
of  the  DNA  in  the  “tails”,  were  not  significantly  different  among  any  of  the  exposure  groups. 
The  conclusion,  that  ML  pulses  were  inducing  multiphoton  absorption  effects  in  the  RPE  cells, 
was  rendered  less  compelling  by  the  inconsistency  of  the  comet  assay  results. 
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Figure  2:  Comet  assays  of  unpigmented  hTERT-RPE  cells  exposed  to  CW  and  ML  810-nm  laser. 
Figure  2A  (left)  shows  results  from  exposures  at  80  W/cm2  (0.5  x  ED50).  Figure  2B  (right)  shows 
results  from  exposures  at  160  W/cm2  (1  x  ED50). 
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DNA  damage  as  a  function  of  pigmentation  load 

Previous  work  in  this  laboratory,1618  as  well  as  by  that  of  Denton  et  al.,11  indicated  that  the 
presence  of  melanin  pigmentation  could  increase  the  susceptibility  of  RPE  cells  to  photo- 
oxidative  stress.  Based  on  these  observations,  experiments  were  performed  to  measure  DNA 
damage  in  cells  with  artificially  enhanced  pigmentation.  RPE  cell  cultures  with  different 
pigmentation  loads  were  produced  by  the  melanosome  supplementation  protocol  described 
above,  and  then  exposed  at  the  ED50  fluence  (160  W/cm2)  to  the  CW  and  ML  laser  protocols. 

Under  these  conditions,  only  cells  with  the  highest  levels  of  melanin  pigmentation,  i.e.  the  cells 
supplemented  with  60  or  180  pi  of  melanosomes,  had  significantly  greater  DNA  damage  as 
assessed  by  comet  tail  length  and  moment  (p<.03,  Bonferroni  multiple  comparison  test).  These 
data  are  shown  in  Figure  3.  Interestingly,  in  the  heavily  pigmented  cells  suffering  the  most  DNA 
strand  breaks,  there  was  no  significant  difference  between  the  comet  assays  in  the  CW  and  ML 
exposed  cells.  This  finding  indicated  that  multiphoton  absorption  was  not  involved  -  or  at  least 
not  exclusively  involved  -  in  the  production  of  the  DNA  strand  breaks  under  these  laser 
exposure  conditions.  These  results  must  be  qualified,  however,  due  to  the  recent  discovery  that 
the  duration  of  the  ML  pulses  was  in  fact  approximately  250  fs.  Multiphoton  absorption 
obviously  becomes  less  likely  for  a  given  irradiance  as  the  pulse  width  increases  and  the  photon 
density  decreases. 


Comet  Assay:  Tail  Length  Comet  Assay.  Tail  Moment 

i&i  cw  EE3  m-l  S52  ew  £Zj  m-l 


Figure  3:  Comet  assay  results  with  HTERT-RPE  cells  pigmented  by  melanosome  supplementation.  All 
exposures  were  made  at  the  ED50  (160  W/cm2).  Figure  3A  (left):  comet  tail  length  measurements.  Figure 
3B  (right):  comet  tail  moment  measurements. 
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In  order  to  emphasize  the  effect  of  cellular  pigmentation  load,  the  results  from  the  melanosome- 
supplemention  experiment  were  replotted  in  Figure  4.  Restricting  the  comet  assay  analysis  to 
tail  length,  the  amount  of  DNA  breakage,  as  a  function  of  the  amount  of  isolated  melanosomes 
added  to  the  cells,  is  shown  in  Figure  4.  Our  previous  work  clearly  showed  that  cellular 
pigmentation  is  directly  related  to  the  amount  of  melanosome  supplementation,  cf.  Figure  2  in 
Glickman  et  al.,  20037.  The  trend  of  the  data  in  Figure  4  indicates  that  there  was  little  or  no 
increase  in  DNA  damage  until  at  least  60  pi  of  melanosomes  were  added  to  the  cells.  While  the 
difference  in  cellular  pigmentation  was  a  major  determinant  in  the  extent  of  DNA  strand 
breakage,  the  laser  output  mode,  i.e.  CW  or  ML,  appeared  to  be  less  important.  This  observation 
motivated  us  to  calculate  the  number  of  photons  absorbed  in  a  typical  RPE  cell  melanosome 
during  the  laser  exposures  used  in  these  experiments,  and,  to  determine  if  sufficient  photons  were 
being  absorbed  to  make  multiphoton  absorption  a  likely  (or  even  possible)  outcome.  These 
calculations  are  developed  in  the  following  sections. 


Effect  of  Melanin  on  DMA  Damage 
(Comet  Assay  Tag  length) 


Figure  4:  DNA  strand  breakage 

following  CW  or  ML  laser,  as  a  function 
of  the  amount  of  cellular  melanin 
pigmentation.  Exposures  were  made  at 
the  EDjo  (160  W/cm2). 


Photon  energy  available  for  absorption  in  melanin 
Initially,  the  energy  available  in  each  photon  is  calculated:19 


E=hi  <» 

This  calculation  yields  2.44  x  10'19  J  as  the  energy  of  each  810  nm  photon.  Next,  the  number  of 
photons  available  in  a  single  laser  pulse  is  calculated  for  the  case  of  a  160  W/cm2  laser  and  a  250 
fs  pulse  duration  (this  was  the  actual  pulsewidth  used  in  the  experiments;  the  planned  pulsewidth 
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of  48  fs  was  not  available).  The  total  energy  per  pulse  is  thus  4  x  10'11  Joules/cm2.  The  number 
of  810  nm  photons  available  per  square  centimeter  is  then  1.6  x  108. 

The  time  for  electron  de-excitation  is  assumed  to  be  less  than  1  x  10'8  seconds  (80  MHz  pulse 
repetition  rate),  so  that  the  energy  from  each  pulse  is  assumed  to  be  non-additive.  Another 
simplifying  assumption  is  that  only  the  physical  area  of  a  melanin  granule  is  available  for 
interaction  with  the  photons.  A  melanin  granule  roughly  approximates  a  rod,  5  microns  in  length 
with  a  1  micron  diameter,20,21  so  the  maximum  surface  area  presenting  for  interaction  would  be 
approximately  a  rectangular  area  of  5  square  microns  (1  micron  x  5  microns  or  5  x  10"8  cm2). 
Thus,  the  maximum  number  of  photons  available  for  primary  interaction  with  a  melanin  granule 
would  be 

5  x  10'8cm2  x  1.6  x  lo8  photons/cm2  =  8  photons  (2) 

This  does  not  account  for  the  cylindrical  nature  of  the  granule,  assumes  all  photons  intersect  the 
granule  along  the  long  axis,  and  should  lead  to  an  overestimate  of  the  number  of  interacting 
photons,  if  area  of  the  granule  is  the  primary  parameter  for  absorption.  Reducing  the  pulse 
irradiance  from  160  W/cm2  or  reducing  the  width  of  the  pulse  will  only  reduce  the  number  of 
photons  available  to  interact  with  the  granule.  A  minimum  of  approximately  200  W/cm2  would 
be  needed  at  48  fs  to  result  in  2  photons  intersecting  a  single  granule. 

Extrapolating  to  the  near  infrared  from  the  visible  spectrum  absorption  coefficient  (pa)  of 
melanin  reported  by  Sardar  et  al.22  yields  a  value  of  0.55  cm'1  for  jia  of  melanin  at  810  nm.  If  a 
simplifying  assumption  is  made  that  the  average  thickness  of  the  cylinder  is  0.5  x  10"4  cm,  then 
using  Beer’s  law,  the  number  of  photons  penetrating  the  granule  would  be 

8(  e-0.55  cm-l  x  0.5  x  10-4  ^  _  y  99973  (3) 

This  would  result  in  2.2  x  10'4  primary  photons  (i.e.  8.0  incident  photons  minus  7.99978  exiting 
photons)  absorbed  per  granule.  This  does  not  account  for  potential  attenuation  of  photons  in 
materials  that  might  result  in  some  energy  transfer  to  the  granules. 
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Photons  required  to  cause  breaks 

Note  that  the  wavelength  for  maximum  biological  effect  of  DNA  in  cells  is  270  nm,  and  3  x  270 
=  810  nm,  implying  that  three  photons  could  possibly  be  effective  at  causing  DNA  damage. 
Assuming  that  the  weak  hydrogen  bonds  associated  with  DNA  require  approximately  3  eV  to 
break,  and  the  energy  of  each  810  nm  photon  is 


E  =  2  44xl0~19 _ ^ _ leV  =  1.525  eV  per  photon. 

photon  1.6x10',9J 


(4) 


In  this  case,  a  minimum  of  2  photons  would  be  required  to  cause  a  break.  If  there  is  100% 
absorption  in  each  granule,  and  it  is  assumed  to  be  absorbed  on  the  side  (rather  than  on  end),  then 
8  photons  would  be  available  to  produce  ionization.  Scaling  this,  only  40  W  would  be  needed  at 
100%  absorption  to  achieve  a  break.  More  realistically,  if  the  absorption  coefficient  for  melanin 
were  used  in  the  calculation,  then  the  minimum  energy  required  for  the  actual  absorption  of  two 
photons  with  total  transfer  of  their  energy  to  the  melanosome  (using  previous  simplifying 
assumptions  regarding  area  and  thickness)  would  be  predicted  as: 

_ 160  W _  2  photons 

2.2x1  O'4  photons  absorbed  5xl0~8cm2  =  2.9  x  1013  W/cm2  (5) 

This  calculation  assumes  that  the  NIR  absorption  coefficient  for  melanin  given  above,  0.55  cm'1, 
applies,  and  only  primary  photons  are  present. 


Range  of  laser  exposures  predicted  to  involve  multiphoton  absorption  as  a  major  damage 
mechanism 

The  number  of  photons  expected  to  impact  a  melanosome  as  a  function  of  laser  radiant  exposure 
can  be  predicted  using  the  methodology  described  in  the  previous  sections.  Some  representative 
calculations  have  been  done  for  a  physiologically  relevant  range  of  exposures,  and  are  presented 
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graphically  below  in  Figures  5  and  6.  Note  that  these  graphs  indicate  only  the  number  of 
photons  incident  on  a  typical,  1  x  5  pm,  RPE  melanosome,  not  the  number  of  photons  actually 
(or  realistically)  absorbed  in  the  pigment  granule. 


Figure  5:  Number  of  photons  incident  on  a  RPE  melanosome  as  a  function  of  irradiance 


Irradiance  (W/cm2) 


Figure  6:  Pulse  duration  required  to  produce  sufficient  irradiance  for  at  least  two  810-nm 
ohotons  to  be  incident  unon  a  RPE  melanosome. 
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DISCUSSION 


Implication  of  results 

Recent  investigation  of  the  principal  damage  mechanisms  induced  in  ocular  tissue  by  ultrashort 
pulse  (<100  fs)  lasers,  have  shown  that  laser-induced  breakdown  has  the  lowest  energy 
threshold,  and  therefore  is  probably  the  predominant  damage  mechanism  at  fluences  near  the 
ED50  for  retinal  and  RPE  lesions. 17,23,24  At  suprathreshold  fluence,  other  mechanisms  are 
activated,  including  photochemical  reactions  producing  oxidative  endpoints.25,26  After  our  initial 
observation  of  increased  DNA  strand  breaks  following  exposure  of  cultured  RPE  cells  to 
ultrashort  pulse  laser,  compared  to  CW  exposures  delivering  the  equivalent  average  power,6  we 
proposed  that  a  damage  mechanism  activated  by  multiphoton  absorption  was  responsible  for  the 
difference.  We  reported  that  two-photon  excitation  of  fluorescence  in  isolated  RPE 
melanosomes  was  due  to  multiphoton  absorption,  but  the  fluorescence  was  observed  at 
suprathreshold  irradiance.5  The  present  investigation  of  DNA  strand  breaks  in  RPE  cells  at 
subthreshold  and  threshold  irradiance  indicated  that  cellular  melanin  plays  a  role  in  mediating 
the  DNA  damage,  but  it  was  not  clear  that  multiphoton  absorption  was  responsible  for  the 
damage  near  the  presumptive  threshold  fluence.  This  conclusion  was  primarily  based  on  the 
finding  of  similar  amounts  of  DNA  damage  produced  by  both  CW  and  ML  laser  exposures  at 
threshold  irradiance.  The  theoretical  analysis  developed  in  this  paper  also  supported  the 
conclusion  that,  at  threshold,  an  insufficient  number  of  photons  are  absorbed  in  a  RPE 
melanosome  to  result  in  multiphoton  absorption. 

Some  qualifications  must  be  considered  in  reaching  this  conclusion.  In  the  analysis  reported 
here,  only  primary  photons,  that  is,  photons  that  travel  directly  from  the  source  to  the  target  are 
considered.  Scattered  photons  are  not  included,  but  are  not  expected  to  contribute  significantly  to 
the  multi-photon  effect.  Another  consideration  is  that  there  may  be  a  specific  “resonance”  effect 
that  causes  a  significant  change  in  the  absorption  properties  of  a  melanin  granule  that  has  yet  to 
be  described.  For  that  matter,  the  optical  properties  of  RPE  melanin  in  the  near  infrared  optical 
radiation  band  have  not  been  precisely  measured.  The  values  used  in  our  analysis  represent  a 
reasonable  extrapolation  from  the  optical  properties  of  melanin  measured  in  the  visible  spectrum. 
Additionally,  the  area  of  interaction  is  assumed  to  be  at  the  most  favorable  presentation  for 
absorption,  and  the  cylindrical  nature  of  the  granule  is  not  included,  so  the  physical  interaction 


14 


area  used  in  these  calculations  is  greater  than  the  actual  physical  size  of  the  granule.  Physical 
size  is  not  always  directly  proportional  to  absorption  coefficients,  and  it  is  known  that  resonance 
reactions  exist  for  many  photon  interactions.19 

Other  mechanisms  besides  direct  photon  absorption  may  also  play  an  important  role  in 
determining  the  damage  to  melanin  granules.  If  melanin  has  undiscovered  or  unusual  properties, 
such  as  the  ability  to  exist  as  a  stable  free  radical,27  then  the  melanosome  may  be  able  to 
maintain  electrons  in  an  excited  state  such  that  high  repetition  rate  photon  pulses  can  produce 
additive  (pulse  to  pulse)  effects.  In  such  a  case,  multiphoton  effects  would  then  be  possible  at 
lower  peak  power  exposures  than  indicated  by  our  analysis.  In  the  absence  of  such  properties, 
however,  individual  pulses  would  require  supra-threshold  energies  on  the  order  of  megawatts  to 
initiate  multi-photon  effects  by  primary  photons. 

Contribution  of  non-linear  interactions  to  possible  multiphoton  effects 

Any  of  the  non-linear  interactions  known  to  occur  with  ultrashort  laser  pulses28'30  could  possibly 
cause  multiphoton  interactions  to  occur.  The  photon  concentrations  produced  could  be  sufficient 
to  cause  localized  ionization  in  each  case.  The  concentration  of  photons  would  be  altered  by  the 
non-linear  reaction,  for  example  through  the  effects  of  self-focusing,  which  might  result  in 
sufficient  energy  deposition  leading  to  alteration  of  absorption  coefficients  with  resulting 
concentration  of  energy.  Such  phenomena  could  manifest  themselves  as  plasma  formation, 
thermal  denaturation  (probably  highly,  localized),  and/or  oxidative  damage. 


Is  multiphoton  absorption  mediated  by  other  possible  cellular  targets? 

Finally,  the  possibility  exists  that  other  intracellular  targets  may  be  involved  in  the  interaction  of 
ultrashort,  infrared  laser  pulses.  In  the  case  of  the  RPE  cell,  we  have  concentrated  on  the 
melansome  because  of  its  prominent  cytoplasmic  location,  its  broadband  optical  absorption 
characteristics,31  and  its  known  involvement  in  ocular  light32  and  laser33  damage.  The  RPE 
melanin  is  clearly  an  important  factor  in  mediating  CW  and  ultrashort  pulse  damage,  as  indicated 
by  the  present  results,  and  it  appears  to  mediate  photo-oxidative  stress  produced  by 
suprathreshold,  ultrashort  pulse  exposures.26  Nevertheless,  at  threshold,  other  cellular  structures 
may  play  a  role  in  cellular  responses  to  ultrashort  laser  pulses.  Mitochondria  are  possible  targets, 
but  -  as  is  also  the  case  with  the  melanosomes  -  it  is  unclear  how,  in  the  short  term,  damage  to 
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these  cytoplasmic  organelles  could  affect  the  nuclear  DNA.  Potential  targets  for  laser  damage 
exist  in  the  nucleus.  As  noted  above,  the  absorption  maximum  for  DNA  damage  occurs  near  270 
nm;  therefore,  a  suitable  absorber  to  mediate  multiphoton  absorption  of  near  infrared  photons  has 
to  be  identified.  The  genomic  DNA  in  eukaryotic  cells  itself  exists  in  nucleosomes,  repeating 
complexes  of  DNA  and  histone  proteins.  The  diameter  of  these  structures  is  on  the  order  of  30 
nm,34  which  is  a  much  smaller  cross-section  compared  to  the  RPE  melanosomes,  and  presumably 
would  absorb  incident  NIR  photons  at  much  lower  efficiency.  The  chain  of  nucleosomes 
comprising  the  DNA,  however,  forms  an  assemblage  that  is  much  larger  than  the  individual 
subunits,  and  if  the  entire  complex  functioned  as  a  single,  absorbing  unit,  then  the  effective 
absorption  cross-section  would  be  correspondingly  larger.  Whether  the  entire  nuclear  DNA- 
histone  complex  functions  as  a  single  absorber  to  promote  multiphoton  absorption  at 
suprathreshold  laser  exposures  remains  to  be  determined.  In  addition,  better  measurements  of 
the  near  infrared  optical  properties  of  cellular  organelles  are  required  in  order  to  evaluate 
hypothetical  non-linear  laser  bioeffects. 
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SECTION  TWO 


Intracellular  signaling  mechanisms  responsive  to  laser-induced 
photochemical  and  thermal  stress 

INTRODUCTION 

The  pigmented  cells  of  the  retinal  pigment  epithelium  (RPE),  by  virtue  of  their  melanin  content, 
absorb  light  over  a  wide  wavelength  range.  With  sufficient  irradiance,  a  variety  of  physiological 
and  physical  responses  including  disruption,  denaturation,  and  photochemical  damage  may 
occur.1'6  The  particular  form  of  the  laser-induced  damage  is  determined  by  the  total  energy 
delivered  to  the  tissue,  the  rate  of  energy  delivery,  the  wavelength  of  the  irradiating  light,  and  the 
optical  properties  of  the  irradiated  tissue.  While  the  response  of  tissue  (especially  ocular  tissue) 
to  laser  exposure  has  been  well  studied  at  the  morphological  level,  the  understanding  of  the 
cellular  response  at  the  molecular  level  is  still  evolving.  For  example,  it  is  becoming  clear  that 
oxidative  stress  may  contribute  to  overall  tissue  damage,  even  if  the  laser  exposure  parameters 
primarily  cause  damage  through  thermal-  or  stress-confined  conditions.7'10  The  present  study 
was  conducted  to  characterize  the  laser-induced  response  of  intracellular  signaling  pathways  that 
are  known  to  be  responsive  to  oxidative  stress.  One  such  signaling  pathway  is  the  transcription 
factor,  NF-kB.11 


The  nuclear  transcription  factor,  NF-kB,  is  a  member  of  a  set  of  proteins  that  controls  the 
expression  of  severaFfamijies  of  genes,  especially  “earlv”  genes  responding  to  stressors,  as  well 
as  those  involved  in  repair  mechanisms.12,13  The  transcription  factor  itself  is  composed  of 
subunits  drawn  from  the  “NF-kB”  class,  designated  p50  and  p52,  as  well  as  from  the  “Rel”  class, 
including  RelA  (also  known  as  p65),  c-Rel,  and  RelB.  These  proteins  share  an  N-terminal 
domain  called  the  Rel  homology  (RH),  that  contains  the  regions  subserving  dimerization  of  the 
protein  subunits,  DNA  binding,  nuclear  localization,  and  binding  of  an  inhibitory  subunit,  IkB.14 
The  most  abundant  form  of  the  transcription  factor  in  the  cytoplasm  of  mammalian  cells  is  a 
dimer  consisting  of  a  light  component,  p50,  and  a  heavy  component,  the  RelA,  or  p65,  subunit. 
This  dimer  is  commonly  referred  to  by  the  generic  name,  NF-kB.15  This  dimer,  when  activated 
by  release  from  its  inhibitory  subunit,  IkB,  is  translocated  into  the  nucleus,  where  it  binds  to 
genomic  DNA  kB  binding  sequences  and  initiates  the  expression  of  select  genes  (Figure  1).  A 
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major  hypothesis  of  the  present  work  is  that  the  melanin  granules  of  the  RPE  cell  represent 
potential  centers  for  interaction  with  incident  light  or  laser  energy,  and  the  absorption  of  this 
energy  results  directly  in  production  of  reactive  species  that  produce  oxidative  damage  to 
cellular  sites.  The  activation  of  NF-kB  appears  to  be  a  response  to  the  laser- induced  photo- 
oxidative  stress. 


stressor  (eg  teser)  ce«»r  targets 


Figure  1.  Schematic  of  NF-kB  regulation,  activation,  degradation,  and  function.  The  figure  illustrates 
the  hypothetical  activation  and  nuclear  translocation  of  NF-kB  following  the  absorption  of  incident  light 
in  the  RPE  melanosomes. 


Regulation  of  NF-kB  activity 

The  basic  regulation  is  that  the  NF-  :;B/Rel  proteins  are  retained  in  the  cytoplasm, 

bound  to  an  inhibitory  subunit,  IkB,  which  blocks  the  nuclear  localization  sequence  of  the  RH 
domain.  This  prevents  the  transcription  factor  from  entering  the  nucleus  until  an  appropriate 
external  signal  or  stressor  causes  dissociation/degradation  of  IkB.  IkB  degradation  is 
accomplished  by  phosphorylation  of  N-terminal  serines  in  IkB  by  a  specialized  kinase,  the  IkB 
kinase  complex  (IKK).  The  second  stage  of  inactivation  involves  ubiquitination  of  IkB  lysine 
sites  through  an  ATP-dependent  reaction.1617  Once  polyubiquitinated,  IkB  undergoes  rapid 
degradation  in  the  proteasome.  Activation  of  the  NF-kB  subunits  can  occur  rapidly,  i.e.,  in 
minutes,  after  a  sufficiently  strong  stimulus  initiates  the  action  of  the  IKK  complex.16  Two  types 
of  IKK  (a  and  P)  have  been  recognized  with  different  and  non-overlapping  functions;  however, 
this  aspect  will  not  be  discussed  further  here.  Following  activation,  the  NF-kB  is  transported 
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into  the  nucleus  by  means  of  the  nuclear  localization  sequence.  Binding  of  the  protein  to  kB 
sequences  in  the  DNA  cause  transcription  of  large  families  of  genes.  Because  there  are  actually 
several  types  of  Rel/NF-KB  protein  subunits  that  can  form  homo-  and  heterodimers,  and  each 
dimer  can  bind  to  different  kB  sites  in  genome,  a  diversity  of  gene  expression  responses  are 
possible,  depending  on  the  particular  dimer  activated.18 

Stressors  causing  activation  of  NF-kB 

Many  external  stressors  activate  NF-kB.  There  is  a  large  literature  describing  NF-kB  activation 
by  toxins,  pathogens  (including  bacteria  and  viruses),  physiological  stress  (e.g.  ischemia  or  shear 
stress),  oxidative  stress,  and  physical  stress  (e.g.  hyperthermia),  as  reviewed  by  Pahl.18  Other 
specific  activators  include  inflammatory  cytokines19,20  and  ionizing  radiation.21’23  With  respect 
to  the  bioeffects  of  light,  there  has  been  much  work  on  the  effect  of  UV  exposure  and  NF-kB 
activation  in  cells,  which  has  found  that  the  generation  of  reactive  oxygen  species  is  a  critical 
precursor  to  activation  of  the  transcription  factor.24'27  Moreover,  UV  exposure  can  activate  NF- 
kB  in  enucleated  cells,  indicating  that  the  sensitive  receptor  element  mediating  the  activation  is 
in  the  cytoplasmic  compartment  -  most  likely  in  the  plasma  membrane.26,28  Changes  in  the 
oxidative  (redox)  state  of  the  cell  are  tightly  linked  to  the  activity  level  of  NF-kB.29  Deliberate 
production  of  reactive  oxygen  species  within  the  cell,  for  example  in  photodynamic  therapy,  is 
an  effective  activator  of  NF-kB.30'32  Consistent  with  the  role  of  oxidative  stress  as  an  activating 
agent,  treatment  with  antioxidants  usually  reduces  the  activation  of  the  transcription  factor.33'35 


Role  of  NF-kB  activation 


The  genes  activated  by  NF-kB  are  involved  in  a  wide  range  of  cellular  activities  and  function. 
About  150  genes  are  known  to  be  activated  by  NF-KB/Rel  proteins.18  These  genes  are  involved 
in  cellular  repair,  restoration  and  proliferation,  as  well  as  production  of  pro-  and  anti-apoptotic 
factors.  Many  of  these  genes  have  contrasting  actions  on  the  cell;  not  all  the  genes  are  activated 
by  every  NF-KB/Rel  family  member.  Particularly  noteworthy  is  that  NF-kB  activation  is 
capable  of  modulating  immune  responses  and  mediating  inflammatory  reactions,  as  evidenced 
by  its  effect  on  the  expression  of  pro-inflammatory  cytokines,  such  as  TNF-a  and  ~p,  IL-la  and 
~P,  IL-6,  and  IFNa.  Recognition  of  the  role  of  NF-kB  in  promoting  or  mediating  an 
inflammatory  response  after  retinal  laser  injuries  would  be  an  important  advance  in 
understanding  the  cellular  responses  to  -  and  possibly  developing  an  effective  therapy  for  -  such 
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injuries.  A  central  question  is  whether  NF-kB  activation  represents  a  death  (apoptotic)  signal,  or 
the  initiation  of  a  repair  response.  The  answer  to  this  question  appears  to  be  complex,  and  highly 
cell  type-dependent.  Depending  on  the  cellular  system,  and  the  type  of  activating  stress 
involved,  NF-kB  may  either  promote  or  prevent  apoptosis.36  In  many  pathological  conditions, 
particularly  in  neurodegenerative  conditions  leading  to  accumulation  of  P-amyloid  precursor  in 
neurons,  the  activation  of  NF-kB  prevents  neuronal  loss  through  apoptosis.37  In  contrast,  NF-kB 
activated  after  ischemia  or  glutamate-induced  neurotoxicity  was  found  to  be  pro-apoptotic.38, 39 
In  ocular  light  damage  studies,  after  photo-oxidative  stress,  retinal  photoreceptors  that  failed  to 
activate  NF-kB  did  not  survive,  undergoing  cell  death  by  either  apoptosis  or  necrosis.40, 41 
Overall,  it  is  probably  more  common  for  NF-kB  to  exert  an  anti-apoptotic  action,  i.e.,  the 
cellular  genetic  response  initiated  by  NF-kB  produces  a  reparative  effect  that  prevents  apoptosis, 
unless  the  external  stressor  exceeds  the  lethal  threshold.36  Regardless  of  the  consequence  of  its 
activation,  NF-kB  serves  as  an  indicator  of  oxidative  stress.  As  noted  above,  its  appearance  in 
the  nucleus  (translocation)  is  a  function  of  oxidative  stress  to  the  cell.  There  is  an  extensive 
literature  on  the  effect  of  generalized  UV  exposure  on  NF-kB  metabolism,  e.g.,  see  review  by 
Legrand-Poels  et  al.42  Characterization  of  NF-kB  activation  following  laser  injury  to  the  retina 
and  retinal  pigment  epithelium  (RPE)  may  be  a  useful  way  to  assess  the  severity  of  the  laser 
injury,  as  well  as  the  efficacy  of  experimental  therapies.  The  research  reported  here  has 
demonstrated  that  NF-kB  is  activated  after  sub- lethal  laser  exposure,  and  the  degree  of  activation 
can  be  modulated  by  certain  antioxidants,  indicating  that  nuclear  NF-kB  levels  are  tied  to  the 
redox  environment  of  the  cell. 


METHODS 
In  vitro  cell  model 

These  experiments  were  performed  on  hTERT-RPE  cells,  a  line  of  human-derived  retinal 
pigment  epithelial  (RPE)  cells,  stably  transfected  with  the  gene  for  human  telomerase  (the 
hTERT-RPE  cell  line  was  obtained  from  Geron  Corporation,  Menlo  Park,  CA).  The  presence  of 
this  gene  prevents  the  onset  of  replicative  senescence  in  these  cells  in  culture;  in  other  respects 
the  cells  exhibit  a  normal  phenotype.43'44  The  cells  were  grown  on  plastic,  24-well  culture  plates, 
and  were  maintained  at  37°  C  in  Dulbecco’s  Minimum  Essential  Medium  (DMEM)  containing 
10%  calf  serum  under  an  atmosphere  of  95%  02  and  5%  C02. 
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Laser  irradiation 

The  hTERT-RPE  cells  were  prepared  for  laser  irradation  by  replacing  the  standard  DMEM 
containing  phenol  red  indicator  with  Dulbecco’s  Phosphate  Buffered  Saline  (DPBS),  which  is 
colorless.  The  culture  plates  were  placed  on  a  computer-controlled  X-Y  stage,  so  that  each  well 
could  be  exposed  to  the  laser  for  a  set  amount  of  time,  or  left  unexposed  as  a  control. 
Continuous  wave  (CW)  exposures  were  made  with  a  Coherent  Model  920  Argon-ion  laser,  the 
output  of  which  consisted  of  a  mixture  of  the  488.1  nm  line  (55%)  and  the  514.5  nm  line  (45%). 
The  culture  plates  were  maintained  at  ambient  room  temperature  during  the  exposure  protocol. 
The  exposure  duration  was  10  m  at  an  irradiance  of  0.5  W/cm2.  Laser  power  was  measured  with 
a  Molectron  PM3  thermopile  detector  driving  an  EP-1000  radiometer.  This  level  of  exposure 
has  previously  been  demonstrated  to  produce  markers  of  photochemcal  stress  in  the  cells,  e.g. 
protein  carbonyl  adducts.45  Following  the  laser  exposure,  the  culture  plates  were  returned  to  the 
37°  C  incubator  until  the  time  points  for  harvesting  the  cells  for  analysis  were  reached. 


NF-kB  Assay 

Cells  were  harvested  at  various  time  points  up  to  24  h  post-exposure.  They  were  freed  from  the 
culture  plate  at  the  desired  time  points  by  replacing  the  DPBS  with  0.25%  trypsin.  After  about  1 
m,  the  cells  were  released  and  transferred  into  plastic  microcentrifuge  tubes.  After  washing 
twice  with  DPBS,  the  cells  were  pelleted  and  frozen  at  -70°  C  until  analysis.  For  NF-kB 
analysis,  nuclear  extracts  from  the  cells  were  subjected  to  an  electrophoretic  mobility  shift  assay 
(EMSA).21  Briefly,  the  cells  were  thavved,  resuspended  iii  lysis  buffer,  and  the  nuclei  were 
separated  from  the  cytoplasmic  lysate  by  centrifugation.  Protein  was  extracted  from  the  nuclei, 
loaded  on  polyacrylamide  gels,  and  probed  with  32P-labeled,  NF-kB -specific  oligonucleotides. 
Binding  to  proteins  in  the  gel  was  imaged  by  autoradiography.  Quantification  of  the  labeled 
bands  in  the  gels  was  accomplished  by  densitometric  analysis  of  the  autoradiographic  images 
using  Image-Pro  software  (Media  Cybernetics,  Silver  Spring,  MD). 


25 


Antioxidant  treatment  of  hTERT-RPE  cells 

Antioxidant  treatment  consisted  of  N-acetylcysteine  (NAC)  at  20  mM,  ascorbic  acid  (AA, 
vitamin  C)  at  500  p.M,  or  pyrrolidine  dithiocarbamate  (PDTC)  at  100  |J.M.  The  antioxidant 
solutions  were  prepared  in  DPBS.  Reagents  were  obtained  from  the  Sigma  Corporation  (St. 
Louis,  MO).  The  concentrations  of  the  antioxidants  were  chosen  based  on  previous 
experimental  observations,  e.g.  with  PDTC  and  NAC  on  human  aortic  endothelial  cells46,  and 
with  AA  on  bovine  RPE  cells.47  The  antioxidants  in  DPBS  were  tested  individually,  in  different 
wells  on  a  culture  plate.  The  cells  were  incubated  with  the  antioxidant  for  30  m  prior  to  laser 
exposure.  Following  this  pretreatment,  the  cells  were  exposed  to  laser,  and  then  prepared  for 
NF-kB  assay,  as  described  above. 


RESULTS 

Induction  of  NF-kB  by  exposure  of  hTERT-RPE  cells  to  CW  visible  laser 
In  non-irradiated  hTERT-RPE  cells,  there  was  very  little  NF-kB  (p50/p65  dimer)  present  (Figure 
2A).  The  p50/p50  dimer  was  present,  and  was  possibly  constitutively  expressed  in  these  cells. 
In  cells  exposed  to  the  CW  argon-ion  laser,  the  p50/p65  dimer  rapidly  appeared  (Figure  2A). 
The  specificity  of  the  assay  was  evaluated  by  incubating  the  nuclear  extracts  with  an  excess  of 


CW  Laser  +  + 
Competitor-  + 


Figure  2.  A  (left):  activation  (nuclear  translocation)  of  p50/p60  dimer  in  RPE  cells  exposed  to  laser  (lane 
2, 60  m  post  exposure  assay),  compared  to  unexposed  controls  (lane  1).  The  amount  of  p50/p50  dimer 
in  the  nucleus  was  not  affected  by  laser  exposure.  B  (right):  Competitive  experiment  showing  specificity 
of  NF-kB  DNA-binding  activity  in  the  EMSA  assay.  Nuclear  extracts  were  incubated  in  the  absence  or 
presence  of  an  excess  of  unlabelled  NF-xB-binding  oligonucleotide  for  5  m,  and  then  probed  with  32P- 
labelled  probe.  Competition  by  the  unlabeled  probe  indicates  specificity  of  assay. 
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unlabeled,  homologous  oligonucleotide  competitor  for  5  m,  running  the  extracts  on 
electrophoretic  gels,  and  then  probing  the  gels  with  32P-labeled  NF-KB-specific  oligonucleotides. 
The  presence  of  the  homologous  competitor  reduced  the  binding  of  the  labelled,  oligonucleotide 
probe  (Figure  2B),  demonstrating  the  specificity  of  the  assay. 

Time-course  of  NF-kB  activation  following  CW  laser  exposure 

CW  argon-ion  laser  exposure  of  hTERT-RPE  cells  induced  a  biphasic  pattern  of  NF-kB 
activation.  After  exposure  to  the  mixed,  488  and  514  nm  output  of  a  CW  Argon-ion  laser  for  10 
m  at  an  irradiance  of  0.5  W/cm2,  an  increase  in  nuclear  NF-kB  was  detected  as  early  as  1  h  post¬ 
exposure,  compared  to  unexposed  cells,  which  exhibited  little  or  no  NF-kB  activation.  The 
activation  followed  a  biphasic  pattern,  with  a  decline  following  the  initial  increase,  followed  by  a 
second,  higher  increase  by  24  h  post-exposure  (Figure  3).  The  early  increase  was  likely  due  to 
activation  of  a  pre-existing  pool  of  cytoplasmic  NF-kB  (i.e.,  released  from  IkB),  while  the  later 
increase  probably  represented  de  novo  synthesis.  Interestingly,  this  biphasic  pattern  was  similar 
to  that  observed  after  exposure  to  ionizing  radiation,21  and  was  also  reported  after  some  types  of 
photosensitized  cell  damage,  with  the  delayed  response  component  beginning  at  about  2  h  post¬ 
treatment  and  continuing  for  24  h.42 


Omin  10  1h  6h  20  h  24  h 

T«mo  after  Exposure 


Figure  3.  Biphasic  activation  of  NF-kB  DNA  binding  activity  indicated  by  EMSA.  RPE  cells  were 
harvested  at  the  indicated  times  after  the  laser  exposure  protocol  described  in  the  text,  and  nuclear  extracts 
were  assayed  for  NF-kB.  The  amount  of  NF-kB  found  at  each  time  point  is  indicated  on  the  graph  as  a 
fold-increase  over  control.  Error  bars  indicate  the  data  range  for  2  experiments. 
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Effects  of  antioxidant  pretreatment  on  laser  activation  of  NF-kB 

Many  reports,  based  on  studies  in  diverse  cell  lines  and  utilizing  a  range  of  external  stressors, 
indicate  that  the  level  of  activation  of  NF-kB  is  modified  by  antioxidant  treatment.48  Two  thiol- 
based  antioxidants  frequently  reported  to  reduce  the  level  of  activated  NF-kB  are  N- 
acetylcysteine  (NAC)  and  pyrrolidine  dithiocarbamate  (PDTC).  The  ability  of  these  agents  to 
inhibit  NF-kB  nuclear  translocation  following  laser  exposure  was  tested  in  the  hTERT-RPE  cell 
line. 

PDTC  was  tested  at  a  concentration  of  100  pM.  The  hTERT-RPE  cells  were  incubated  at  37°  C 
for  30  m  in  DPBS  containing  100  pM  PDTC,  and  then  exposed  to  CW  argon-ion  laser  as 
described  above.  Cells  were  harvested  for  NF-kB  assay  at  1  h  and  24  h  post  laser  exposure. 
Pretreatment  with  PDTC  failed  to  reduce  the  activation  of  NF-kB;  indeed,  the  level  of  the  factor 
in  nuclear  extracts  was  increased  after  PDTC,  even  in  the  absence  of  laser  exposure  (Figure  4). 


Figure  4.  PDTC  pretreatment  of  RPE  cells  does  not  reduce 
NF-kB  activation  by  laser.  RPE  cells  were  incubated  in  1 00 
pM  PDTC  prior  to  laser  exposure.  NF-kB  was  assayed  at 
24  h  and  no  change  was  found  in  the  amount  of  p50/p65 
dimer  activated  by  the  laser,  compared  to  untreated  cells. 
In  addition,  unexposed  cells  treated  with  PDTC  exhibited 
higher  levels  of  nuclear  p50/p65  dimer.  The  amount  of 
p50/p50  dimer  in  laser-exposed  cells  was  slightly  reduced 
by  the  PDTC  treatment.  Results  of  densitometxic  scans  are 
shown  below  each  lane  in  the  EMSA  assay  gel.  Upper  bars: 
p50/p65;  lower  bars:  p50/p50. 


Similarly,  pretreatment  of  the  hTERT-RPE  cells  with  NAC  (20  mM  in  DPBS)  did  not  reduce  the 
laser-induced  activation  of  NF-kB.  At  the  1  h  and  the  24  h  assay  time  points,  the  laser-exposed 
cells  treated  with  NAC  had  higher  levels  of  nuclear  NF-kB,  compared  to  the  untreated  cells 
(Figure  5).  Thus,  the  two  thiol-based  antioxidants  did  not  reduce  the  stimulatory  effect  of  the 
laser  on  the  NF-kB  signaling  pathway,  and  possibly  enhanced  it. 
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Figure  5.  NAC  pretreatment  of  RPE  cells  does  not 
reduce  NF-kB  activation  by  laser.  RPE  cells  were 
incubated  in  20  mM  NAC  prior  to  laser  exposure.  This 
pretreatment  did  not  reduce  NF-kB  activation  either  at 
1  h  or  at  24  hr.  In  the  absence  of  laser  exposure,  the 
NAC  treatment  increased  the  amount  of  NF-kB  in  the 
nuclear  extracts.  Results  of  densitometric  scans  are 
shown  below  each  lane  in  the  EMSA  gel  (upper  bars: 


In  contrast  to  the  two  thiol-based  antioxidants,  the  30  m  pretreatment  of  the  hTERT-RPE  cells 
with  AA  (500  pM  in  DPBS)  reduced  the  activation  of  NF-kB  by  laser  exposure.  This  reduction 
by  AA  was  evident  at  the  6  h  and  24  h  post-exposure  assay  time  points  (Figure  6).  At  6  h,  the 
reduction  in  NF-kB  was  43%,  and  at  24  h,  the  reduction  was  75%,  compared  to  untreated,  laser- 
exposed  cells.  Interestingly,  in  cells  not  exposed  to  laser,  but  treated  with  AA,  there  was  a  slight 
increase  in  the  level  of  activated  NF-kB:  approximately  a  20%  increase  at  6  h,  and  a  3-fold 
increase  at  24  h. 
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Figure  6.  AA  pretreatment  of  RPE  cells  reduces  NF-kB 
activation  by  laser.  RPE  cells  were  incubated  in  500  pM 
AA  prior  to  laser  exposure.  The  amount  of  NF-kB  in  the 
nuclear  extracts  was  reduced,  especially  at  the  24  h  time 
point.  AA  treatment  by  itself  appeared  to  increase  the 
amount  of  p50/p65  dimer  in  the  control  cells.  Results  of 
densitometric  scans  are  shown  below  each  lane  in  the 
EMSA  gel  (upper  set  of  bars:  p50/p65  dimer;  lower  set  of 
bars:  p50/p50  dimer). 
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DISCUSSION 


NF-kB  as  an  indicator  of  photo-oxidative  and  thermal  stress 

The  CW  laser  irradiation  used  in  this  study  stressed  the  hTERT-RPE  cells,  as  evidenced  by  the 
activation  of  NF-kB,  which  is  generally  an  indicator  of  oxidative  stress  or  a  response  to  a  pro- 
inflammatory  signal,  such  as  a  cytokine.48  The  laser  irradiance  used  in  these  experiments  has 
been  demonstrated  previously  to  produce  endpoints  associated  with  oxidative  damage,  including 
depletion  of  antioxidants,49  carbonyl  adducts  to  cellular  proteins,45  and  peroxidation  of  fatty 
acids.50  Although  the  laser  exposure  paradigm  was  designed  to  minimize  temperature  transients 
in  the  target  RPE  cells,  the  question  of  the  effect  of  thermal  stress  on  the  cells  may  be  raised.  In 
general,  hyperthermia,  or  heat  shock,  has  been  found  to  decrease  or  inhibit  the  activation  of  NF- 
kB,51,  52  although  some  neural  tissues  have  been  reported  to  show  an  increase.53  Heat  shock 
may  reduce  the  levels  of  NF-kB  by  producing  insolubilization  of  IkK.54  Nevertheless,  because 
the  moderate  laser  exposure  irradiance  used  in  this  work,  0.5  W/cm2,  was  below  the  level 
generally  considered  to  produce  thermal  tissue  damage,  and  yet  produced  an  upregulation  of  NF- 
kB  in  RPE  cells,  we  concluded  that  the  primary  effect  of  the  laser  exposure  was  to  induce  photo- 
oxidative  stress. 

Lack  of  effect  of  thiol-based  antioxidants 

A  surprising  finding  was  that  the  thiol  antioxidants,  NAC  and  PDTC,  did  not  inhibit  the  NF-kB 
activation  caused  by  laser  exposure.  Many  reports  have  found  that  one  or  both  of  these  agents, 
as  well  as  other  thioreducing  agents,  interfere  with  part  of  the  transcription  factor  signaling 
cascade.34,55  The  ability  of  antioxidants  to  prevent  activation  of  NF-kB  has  led  to  the  general 
hypothesis  that  reactive  oxygen  species  (ROS)  are  somehow  involved  in  triggering  the 
phosphorylation  and  degradation  of  the  IkB  inhibitory  subunit,  thereby  releasing  the  active 
dimer  component.29  It  is  to  be  noted,  however,  in  many  studies  the  inducing  agent  for  NF-kB  is 
a  pro-inflammatory  cytokine,  such  as  tumor  necrosis  factor  (TNFa),  an  oxidizing  chemical,  or 
other  agent  that  directly  interacts  with  receptor  proteins  in  the  plasma  cell  membrane.56 
Although  thiol  reducing  agents  were  supposed  to  have  prevented  the  generation  of  ROS,  a  recent 
study  has  shown  that  the  likely  mechanisms  of  action  of  these  agents  are  unrelated  to  their 
antioxidant  activity.  For  example,  NAC  was  shown  to  reduce  the  binding  efficiency  of  TNFa 
with  its  membrane  receptor,  and  PDTC  reduced  the  efficiency  of  ubiquitin  ligase,  which  is  an 
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enzyme  involved  in  the  inactivation  of  IkB.57  The  lack  of  effect  of  the  thiol  reducing  agents 
against  photo-oxidative  stress  indicates  that  the  laser  radiation  is  interacting  with  yet  another 
cellular  component  to  produce  the  transcription  factor  response. 

The  RPE  cell  melanosome  as  the  site  for  generation  of  photo-oxidative  stress 
In  contrast  to  the  thiol-based  antioxidants,  AA  pretreatment  was  found  to  be  highly  effective  in 
reducing  activation  of  NF-kB.  This  observation,  combined  with  the  fact  that  the  hTERT-RPE 
cells  are  able  to  synthesize  melanin  pigment  in  culture,58  raises  the  possibility  that  the  pigment 
granules  are  the  site  of  oxidative  stress  in  these  cells.  At  the  time  these  experiments  were 
conducted,  the  hTERT-RPE  cells  were  noted  to  have  pigment  granules  (pigmentation  was  visible 
in  the  cell  pellets).  The  ability  of  AA,  but  not  thiol-based  antioxidants,  to  quench  light-activated 
free  radical  reactions  of  melanin  was  noted  in  previous  work.47  The  failure  of  thiol-based 
antioxidants  was  ascribed  to  the  higher  oxidation  potential  of  the  -SH  (thiol)  group,  +920  mV, 
compared  to  that  of  ascorbic  acid,  +320  mV.  Thus,  if  the  action  of  the  Argon-ion  CW  laser 
exposure  was  to  excite  reactive  melanin  free  radicals  which,  left  unquenched,  would  produce 
oxidative  damage  at  specific  cellular  sites,  then  AA  would  be  expected  to  reduce  these  effects, 
while  the  thiol-based  antioxidants  would  not.  These  observations  support  our  hypothesis  that  the 
cytoplasmic  pigment  granules  are  the  locus  for  laser-induced  oxidative  interactions  in  the  RPE 
cell. 


The  ability  of  AA  to  reduce  the  activation  of  NF-kB  also  raises  an  interesting  speculation  about 
•the  difference  between  oxidative  stress  produced  by  molecular  or  chemical  agents,  and  that 
produced  by  light  exposure.  In  the  case  of  TNFcc -mediated  induction  of  NF-kB,  the  presumptive 
ROS  are  produced  by  the  stimulation  of  biochemical  pathways,  such  as  Rac/NADPH  oxidase, 
Rac  being  a  GTP-binding  protein  involved  in  many  signaling  cascades.57  The  production  of 
ROS  after  TNFa  involves  activation  of  energy-dependent  enzymatic  reactions.  In  contrast,  the 
direct  effect  of  light  on  melanin  to  produce  reactive  free  radicals  does  not  require  additional 
energy,  nor  are  other  enzymatic  pathways  necessarily  involved  (but  may  be  affected).50  Thus, 
suitable  antioxidants  may  directly  quench  ROS  produced  by  photo-oxidative  mechanisms,  and 
may  potentially  have  a  therapeutic  role  in  reducing  light-induced  damage  in  tissues,  such  as  the 
retina  and  RPE.  It  is  yet  unclear  if  NF-kB  plays  a  pro-  or  anti-apoptotic  role  in  the  RPE.  Recent 
work  has  shown  that  RPE  cells  (not  hTERT-RPE,  but  another  human-derived  cell  line)  remain 
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resistant  to  TNFa-induced  apoptosis,  despite  carrying  a  mutant,  non-degrading  IkB  that  strongly 
inhibits  NF-kB.59  In  the  present  work,  the  hTERT-RPE  cells  exhibited  very  little  apoptosis 
following  the  laser  exposure.  Whether  these  cells’  survival  was  due  to  activation  of  NF-kB,  or 
to  some  other  factor,  remains  to  be  demonstrated. 
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SECTION  THREE 


Non-invasive  thermography  of  laser-induced  hyperthermia  using  magnetic  resonance 

INTRODUCTION 

Recently,  there  has  been  renewed  interest  in  using  laser  heating  to  induce  localized  hyperthermia  to 
treat  cancer  and  other  diseases.1'3  Local  hyperthermia  can  be  used  for  a  wide  variety  of  medical 
interventions  such  as  heart  arrhythmias4  and  treatment  of  solid  and  vascular  tumors,  as  well  as  other 
vascular  anomalies  of  various  etiology  of  the  retina  and  choroid.5  Generally,  photothermal  therapy 
does  not  produce  as  much  widespread  damage  as  does  radiation  therapy,  e.g.  gamma  and  beta 
radiation.6  Therefore  photothermal  therapy  has  been  recommended  as  the  treatment  method  of 
choice  for  small  melanomas  situated  close  to  the  macula  or  optic  nerve.5  Such  treatment  of  tumors 
in  the  eye  (e.g.  retinoblastomas  and  malignant  melanomas),  utilizing  near  infrared  laser  irradiation 
at  the  target  tissue  through  the  pupil,  is  called  Transpupillary  Thermotherapy  (TTT),  which  was 
introduced  by  Oosterhuis.7  The  near  infrared  range  of  electromagnetic  radiation  800-810  nm  laser, 
is  the  best  for  such  procedures  because  of  the  good  optical  penetration  properties  of  this  radiation  in 
the  ocular  tissue.6 

The  object  of  thermal  treatment  is  to  provide  adequate  heating  to  destroy  the  entire  target  volume 
while  sparing  adjacent  normal  structures.  Interestingly,  it  has  been  reported  that  tumor  tissue  is 
more  sensitive  to  increased  temperature  than  is  normal  tissue.8  In  the  TTT  treatment,  strict 
.bonfiferbafr;  of  the  extent  of  photothefmaj  .dj$nSgeisv  entibafrbeeause  lesions  are  frequently 
located  close  to  the  macula  or  optic  nerve  as  well  as  retinal  pigment  epithelium  and  choroid.6  Heat 
conduction  through  diffusion  and  perfusion  processes  may  vary  locally  as  a  function  of  tissue 
architecture,  tissue  composition,  local  physiological  parameters,  and  temperature  itself. 

Thus,  there  is  a  need  for  optimal  photothermal  treatment  where  optimal  irradiation  therapy  should 
take  into  account  the  distribution  of  both  radiant  and  thermal  energy  throughout  the  volume  of  the 
treated  organ  or  tissue.  Optimizing  the  variable  parameters  will  maximize  the  effectiveness  of  the 
hyeperthermia  treatment.  Magnetic  resonance  imaging  (MRI)  has  recently  enabled  a  noninvasive 
method  to  measure  temperature  changes  in  irradiated  tissue,  which  is  more  suitable  and 
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advantageous  in  clinical  use  than  conventional  feedback  temperature  monitoring  using  a  series  of 
thermistor  probes.  It  also  minimizes  stress  on  the  patient,  and  reduces  the  possibility  of 
confounding  interactions  between  the  laser  radiation  and  the  temperature  sensors.9'11  Magnetic 
resonance  thermal  imaging  depends  on  several  parameters  that  are  known  to  be  temperature 
dependent  such  as  the  relaxation  time,  the  diffusion  coefficient,  and  the  proton  resonance 
frequency.12,13  The  excellent  linearity  and  near-independence  with  respect  to  tissue  type,  together 
with  good  temperature  sensitivity,  make  proton  resonance  frequency  (PRF)  based  temperature  MRI 
the  preferred  choice  of  many  applications. 

In  this  investigation  we  carried  out  measurements  of  laser-induced  temperature  changes  in  an  ocular 
phantom  using  the  PRF  based  method  to  evaluate  its  suitability  for  ophthalmic  applications.  This 
method  was  able  to  produce  three-dimensional  temperature  measurements  with  sufficient 
spatiotemporal  resolution  for  practical  measurements  in  the  eye.  Because  of  local  differences  in  heat 
conduction  and  energy  absorption,  three-dimensional  real-time  temperature  mapping  would  be 
ideal.  A  numerical  model  based  on  classical  solution  of  the  heat  diffusion  equation  was  used  to 
calculate  temperature  distribution  in  the  phantom  irradiated  by  a  laser  beam  at  near  infrared 
wavelength  (806  nm),  taking  the  redistributive  scattering  into  account,  which  was  then  compared  to 
measured  empirical  data. 


METHODS 

— — - — •. 

Two  types  of  phantoms  were  used  during  the  course  of  this  research.  A  uniform  ocular  phantom 
was  constructed  by  filling  a  22-mm  diameter,  cylindrical  glass  cuvette  with  either  clear,  1%  agarose 
(Sigma  Chemical  Co.,  St.  Louis,  MO),  or  1%  agarose  containing  a  uniform  suspension  of 
melanosomes  isolated  from  bovine  retinal  pigment  epithelium  (RPE)  as  previously  described.14  The 
agarose  was  dissolved  in  distilled  water,  heated  to  60  °C,  poured  into  the  cuvette,  and  allowed  to 
gel.  The  uniform  phantoms  were  used  only  for  calibration  of  the  MRT  temperature  response,  as 
described  below  (“MR  Thermography”).  For  experiments  involving  laser  radiation,  a  non-uniform 
phantom  containing  a  thin  layer  of  melanin  was  constructed  using  a  10-cc  Pyrex  glass  beaker,  also 
of  approximately  22-mm  diameter.  The  melanin  layer  was  formed  by  mixing  a  suspension  of 
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bovine  RPE  melanosomes  in  agarose  at  a  granule  density  estimated  to  be  similar  to  that  in  native 
RPE,  i.e.  approximately  1010  melanin  gran/ml.14  This  melanosome-agarose  suspension  was 
sandwiched  and  allowed  to  gel  between  two  microscope  cover  slips,  forming  a  layer  of  1 .5  mm 
thickness.  The  sandwich  was  placed  in  the  middle  of  the  beaker,  and  clear  agarose  was  poured  in 
around  it. 


Infrared  Laser  Irradiation  of  Phantoms 

The  laser  source  used  in  this  study  was  an  HPM-F1  laser  diode  module  with  a  built-in  fiber  launch, 
that  produced  1  watt  CW  maximum  output  power  at  806  nm  (Power  Technology  Corp,  Little  Rock, 
AK).  The  experimental  optical  and  temperature-recording  arrangement  is  shown  in  Figure  1.  The 
laser  console  was  situated  about  4.5  meters  away  from  the  center  of  the  MR  magnetic  field.  The 
laser  emission  was  delivered  into  the  cavity  of  the  MR  instrument  using  a  5-meter  long,  600  pm 
diameter  core,  multimode  optical  fiber  (Oz  Optics,  Carp,  Ontario,  Canada).  An  in-line,  fiber  optic 
variable  attenuator  (Oz  Optics)  controlled  the  laser  power  delivered  to  the  phantom.  Two  lenses 
were  used  to  collimate  (LI,  50  mm  focal  length)  and  focus  (L2,  1 1  mm  focal  length)  the  laser  beam 
on  the  phantom  target  inside  the  MR  cavity,  as  shown  in  Figure  1.  The  lenses  and  fiber  optic  were 
mounted  and  aligned  in  a  non-magnetic  holder.  Separately,  the  focused  laser  spot  size  at  the  target 
was  imaged  with  a  BeamView  Analyzer  PC  laser  beam  profiler  (Coherent  Inc.,  Instruments  Div., 
Auburn,  CA  )  and  was  found  to  be  an  approximately  Gaussian  beam.  In  these  experiments,  two 
-laser  bear"  sizes  on  the  target  were  used:  a  small  spot  with  I/s2  radius  of  0.8  mm,  and  a  large  spot 
with  1/e2  radius  of  2.4  mm. 

The  laser  power  at  the  target  position  was  measured  with  an  EPM  1000  energy/power  meter 
equipped  with  a  MP3  thermopile  detector  head  (Molectron,  Eugene  OR).  Depending  on  the 
experiment,  the  target  irradiance  was  varied  over  the  range  of  2.4  -  21.6  W/cm2.  This  was 
determined  by  placing  the  power  detector  at  the  position  of  the  phantom,  and  taking  in  account  the 
small  attenuation  due  to  the  glass  vessel  and  the  agarose.  Laser  beam  alignment  was  visually 
checked  using  a  Find-R-Scope  infrared  viewer  (FJW  Optical  Systems,  Palatine,  IL)  to  ensure  that 
the  laser  beam  was  centered  in  the  ocular  phantom. 
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MR  Thermography  _  . .  . . 

All  MRT  experiments  were  performed  in  the  MRI  Laboratory  of  the  Research  Imaging  Center  at  the 
University  of  Texas  Health  Science  Center  at  San  Antonio.  Measurements  were  made  using  a  PRF 
method  similar  to  that  described  by  Vitkin  et  al.15  For  MR  thermography  methods  utilizing  PRF, 
the  image  phase  difference  is  linearly  dependent  on  the  change  in  temperature  as  follows: 

A(p  =  a  •  y  •  Bo  •  TE  •  AT  (1) 

where  Aq> :  phase  difference,  y  :  the  gyromagnetic  ratio  for  protons,  Bo  :  magnetic  field  strength,  a  : 
the  thermal  coefficient  (i.e.  the  chemical  shift),  TE:  echo  time,  and  AT  :  change  in  temperature.  A 
two-turn  radioffequency  surface  coil  of  6-cm  diameter  was  used  as  the  transceiver.  Images  were 
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obtained  using  an  Oxford  Magnet  Technology  2T,  45-cm  diameter  bore,  superconducting  magnet 
(Oxfordshire,  England),  equipped  with  actively  shielded  Magnex  gradients  and  a  Tecmag  Libra 
console  (Tecmag,  Houston,  TX). 

In  order  to  calibrate  the  MRT  measurements  in  terms  of  absolute  temperature,  an  experiment  was 
performed  with  a  phantom  containing  a  uniform  distribution  of  melanin  granules,  constructed  as 
described  above  (“Ocular  Phantoms”).  The  phantom  was  heated  by  passing  warm  water  from  an 
external  circulating  water  bath  through  Tyvek  tubing  wrapped  around  the  phantom.  When  thermal 
equilibrium  was  reached,  the  water  flow  was  stopped.  MRT  was  performed  by  subtracting  the 
phase  images  of  the  ocular  phantom  at  two  different  temperatures,  and  constructing  a  phase 
difference  map,  from  which  the  temperature  difference  was  calculated  using  equation  (1).  To 
correct  for  phase  drifts  due  to  ambient  temperature  changes  unrelated  to  temperature  shifts  in  the 
ocular  phantom,  thermally  isolated  reference  phantoms,  consisting  of  small  tubes  of  agarose,  were 
placed  adjacent  to  the  ocular  phantom.  On  the  assumption  that  the  phase  in  the  image  of  the 
reference  phantoms  remained  constant,  the  relative  temperature  change  map  of  the  ocular  phantom 
was  corrected  by  comparison  to  the  phase  of  three  points  in  the  image  of  the  reference  phantoms 
(three  points  allowed  a  image  plane  to  be  defined).  The  reference  phantoms  had  known  length, 
which  enabled  the  correction  of  three-dimensional  images.  Gradient  echo  images  were  obtained  at 
various  phantom  temperatures.  From  these  images,  the  value  of  a  (the  chemical  shift)  was  found  to 
be  0.0102  ppm/°C,  consistent  with  general  PRF  temperature  measurements  for  biological  tissues.12 
The  signal  to  noise  ratio  (SNR)  for  the  temperature  measurements  was  calculated  as  the  inverse  of 
the  standard  deviation  of  the  phase  over  the  region  of  interest  (ROI)  in  the  heated  sample.  For 
absolute  calibration,  the  temperatures  inside  the  ocular  phantom  and  in  the  references  were 
measured  independently,  with  an  accuracy  of  0.1  °C,  using  Vitek  fiber  optic  temperature  probes 
(Boulder,  CO)  placed  in  the  phantoms. 

For  all  laser-induced  heating  experiments,  the  layered  phantom  consisting  of  a  melanin  layer 
embedded  in  clear  agarose  as  described  above,  was  used.  The  melanin  itself  has  natural  T1  contrast 
in  MRI.  The  806  nm  light  emitted  from  the  diode  laser  source  was  not  significantly  absorbed  in  the 
clear  agarose,  but  was  almost  entirely  absorbed  and  scattered  within  the  melanosome  layer  (refer  to 
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Figure  3A).  This  was  considered  to  be  a  close  approximation  of  the  behavior  of  near  infrared  light 
in  the  vertebrate  eye.  To  obtain  a  time  course  of  the  heating  in  the  phantom,  a  two-dimensional, 
gradient  echo  imaging  sequence  was  used  with  repetition  time  (TR)  of  210  ms,  echo  time  (TE)  of 
12  ms,  and  the  flip  angle  was  40°.  Eight  slices  were  acquired  with  images  obtained  every  22.9  s. 
The  in-plane  resolution  was  0.25  mm  with  a  slice  thickness  of  1.5  mm.  Six  images  were  taken 
before  the  laser  was  turned  on  and  these  were  averaged  to  provide  a  phase  reference  image.  Laser- 
induced  heating  was  applied  for  15  -  20  minutes,  and  images  were  taken  sequentially  throughout  the 
exposure  period. 

A  three-dimensional  (3D)  temperature  map  was  obtained  in  the  melanin-layer  phantom,  using  a  3D 
gradient  echo  imaging  sequence  providing  an  image  of  0.25x0.25x0.25  mm  resolution  in  each 
direction.  In  this  gradient  sequence,  TR  =  40  ms,  TE  =  12  ms,  and  the  flip  angle  =  20°.  The  time 
required  for  a  total  3D  scan  was  8  minutes.  A  3D  image  prior  to  heating  was  acquired  as  the 
reference  image,  as  well  as  two  3D  images  after  the  laser  had  been  turned  on.  Acquisition  of  the 
first  image  commenced  with  the  onset  of  laser  exposure  (the  center  of  k-space  was  after  4  minutes), 
while  the  second  image  was  acquired  after  the  completion  of  the  laser  exposure  (the  center  of  k- 
space  was  after  13  minutes).  Three-dimensional  scans  were  acquired  during  heating  of  the  phantom 
with  the  small  (0.8  mm  radius)  as  well  as  the  large  (2.4  mm  radius)  laser  beam  spot  sizes. 

Image  Analysis 

Phase  maps  were  reconstructed  from  adjusted  phase  drifts  using  IDL  (Interactive  Data  Language 
from  Research  Systems,  Inc.,  Boulder,  CO),  using  the  three-point  phase  correction  method 
described  in  section  2.3,  and  by  previous  authors.11,22  Temperature  maps  were  then  calculated  using 
equation  (1). 

Infrared  Imaging 

Infrared  photography  was  performed  using  a  Sony  F717  digital  camera  set  to  its  “Nightshooting” 
operational  mode,  i.e.  with  the  CCD  imager  set  to  record  only  near  infrared  light.  Images  were 
taken  of  the  phantom  during  806  nm  laser  irradiation,  and  then  transferred  to  a  personal  computer 
for  processing  and  display. 
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Thermal  Model 


A  basic  model  for  the  transient  thermal  tissue  response  to  continuous  diode  laser  output  was  used  to 
predict  the  temperature  rise  and  to  compare  it  with  the  MRT  measured  values.  The  thermal  model 
accounted  for  the  effect  of  light  scattering,  which  is  predominant  in  the  case  of  melanized  tissue.  In 
previous  work,  the  inverse  adding  doubling  (IAD)  method16  was  used  to  determine  optical 
properties  of  RPE  melanin.  The  scattering  coefficients  were  found  to  be  much  higher  than  the 
absorption  coefficients  of  melanin  at  visible  wavelengths,17,18  and  this  relationship  would  be 
expected  to  be  similar  in  the  near  infrared,  although  the  relative  magnitude  of  each  coefficient 
would  be  lower  than  the  corresponding  visible  wavelength  values.18  Moreover,  the  exact  values 
depend  slightly  on  the  concentration  of  melanin  in  the  phantom.  Thus,  when  compared  with  visible 
wavelength  lasers,  lasers  emitting  in  the  near  infrared  (e.g.  806  nm  wavelength)  will  experience  less 
scattering  and  correspondingly  increased  transmission  through  non-transparent  (turbid)  media. 

Some  simplifying  assumptions  were  made.  The  thermal  radiation  away  from  the  phantom  walls 
was  neglected,  and  material  constants  were  taken  to  be  uniform  and  time  independent.  The  optical 
and  properties  of  the  material  were  also  considered  to  remain  constant  during  temperature  elevation. 

If  the  incident  intensity,  I0(r)  (Watt/cm2),  of  a  cylindrical,  symmetric  laser  beam  penetrates  the 
surface,  then  the  intensity  of  the  direct  beam,  IB  follows  as  in  Beer’s  law: 

I®(r,z)  =  I0(r)exp(-//z)  . .  _  . . . .  . . . . (2) 

where 

M  =  Ma+Ms  (3) 

and  p  is  the  extinction  coefficient,  pa  is  the  (direct)  absorption  coefficient  of  radiation,  and  ps  is 
scattering  coefficient  (redistribution)  of  radiation  out  of  the  beam.  However,  the  intensity  of 
scattering  radiation,  Is,  at  (r,z)  is  the  sum  of  direct  intensity  as  result  of  depolarization  and  other 
randomizing  effects,  and  can  be  written  as  for  cylindrical  symmetry19,20: 
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(4) 


s(r,z)  =  Jr'dr'  jd0"  jdz'  x 


I0(r')exP(-Mz'-lMs) 
l 2 


where  /  is  the  distance  between  points  (i.e.  intensity  of  radiation  arriving  at  (r',0',z'),  to  which  light 
has  been  scattered  from  (r,0,z)  with  coefficient  ps  in  cylindrical  co-ordinates). 

From  the  input  optical  parameters,  the  effective  penetration  depth  5eff  can  be  written: 


where  g  is  the  anisotropy  factor  (dimensionless). 

The  temperature  response  of  tissue  to  a  laser  radiation  source  S  is  governed  by  the  following  heat 
equation:21 


pc 


—  =  KV2T+S 
at 


(6) 


where  T  is  temperature,  p  is  density  (kg/m3),  c  is  the  specific  heat  (J/kg  •  °K)  of  the  substance,  and 
K  is  the  thermal  conductivity  (W/m  •  °K).  The  source  function  can  be  expressed  as: 

S(r,z,t)  =  a[IB(r,z,t)  +  r(r,z,t)]  (7) 

This  enters  the  heat  equation  linearly  as  a  driving  function.  The  laser  beam  profile  is  considered  as 
a  Gaussian  profile  with  laser  beam  waist  r0: 


I0(r)  =  h  exp[-(r  I ra)2~\ 


(8) 


This  distribution  corresponds  to  total  input  power  : 


P  =  xr02Io  (9) 

An  analytical  solution  of  the  heat  differential  equation  by  using  the  Eigenfunction  with  the  relation 
between  /  and  f  has  been  discussed  in  detail.19  In  the  present  application,  these  problems  were 
solved  numerically  by  using  Mathcad  (Mathsoft  Engineering,  Cambridge,  MA);  the  results  of  the 
computation  were  displayed  graphically,  as  will  be  shown  below. 


RESULTS 

Measurement  of  Heating  Time  course 

Temperature  profiles  within  the  layered  phantom  during  continuous  application  of  the  806  nm  diode 
laser  were  measured  using  MRT,  and  compared  to  the  theoretical  values  predicted  from  the  heat 
diffusion  model.  The  optical  properties  assigned  to  the  melanin  layer  phantom  were  as  follows: 
(scattering  coefficient)  ps  -11.4  cm"1,  (absorption  coefficient)  pa  -0.55  cm"1,  and  (anisotropy  factor) 
g  =  0.92.  The  thermal  constants  were  taken  as  pc=  4.4  J/cm3  •  °C  (i.e.  as  for  water),  and  K=7xl0'3 
W/cm  •  °C,  and  were  assumed  to  be  uniform  throughout  the  phantom. 

During  steady  exposure,  the  time  course  of  heating  was  obtained  by  serial  MRT  measurements  as 
described  above  (“MR  .Thermography”)  ..The  layered  phantom  reached  15%  thermal  equilibrium, 
after  a  laser  exposure  of  25-30  seconds,  as  shown  in  Figure  2,  with  a  target  irradiance  of  21.6 
W/cm2  at  the  center  of  the  phantom,  and  the  0.8  mm  radius  (1/e2)  laser  beam.  The  experimental 
example  shows  that  the  steady  state  temperature  reached  within  75  sec  was  23  °C  above  the  initial 
baseline,  which  is  very  high  relative  to  the  temperature  achieved  in  typical  hyperthermia  treatment. 
Also  shown  in  Figure  2  is  an  estimate  of  the  error  in  the  two-dimensional  MRT  measurements. 
Prior  to  the  start  of  laser  heating,  the  standard  deviation  of  the  phase  in  six  consecutive  2D  images 
was  0.047  radians,  corresponding  to  0.72  °C.  This  implied  that  the  temperature  resolution  in  the 
unheated  sample  was  better  than  1  degree.  During  laser  heating,  however,  the  signal  to  noise  ratio 
declined,  and  the  standard  deviation  of  the  phase  increased,  leading  to  an  estimate  of  measurement 
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error  of  1.4  °C.  This  is  indicated  by  the  error  bars  in  Figure  2  for  the  time  points  starting  at  75  sec 
after  the  onset  of  laser  exposure,  i.e.  for  MRT  measurements  made  during  steady-state  heating. 


Figure  2.  Heating  time  course  in  the  layered  phantom  during  CW  infrared  laser  exposure,  as 
measured  with  MRT.  The  solid  line  shows  the  temperature  change  at  the  hottest  point  in  the  phantom, 
within  the  melanin  layer,  measured  with  MRT  in  sequential  images  during  laser  irradiation.  The 
dotted  line  shows  theheating  time  course  predicted  by  the  heat  diffusion  model.  The  error  bars  on  the 
MRT  measurements  show  the  standard  deviation  of  the  phase  measurements,  converted  to 
temperatures.  Laser  heating  was  applied  for  15  min;  the  temperature  measured  in  the  phantom  during 
the  first  400  sec  of  the  experiment  is  shown  in  this  figure.  The  ordinate  shows  the  temperature 
increase  in  the  phantom’s  melanin  layer  during  laser  irradiation.  Taking  into  account  minor  scattering 
and  absorption  of  the  806  nm  laser  in  the  clear  agarose  portion  of  the  phantom,  the  target  irradiance  in 
the  melanin  laver  was  calculated  as  21.6  W.cm2.  where  the  laser  beam  radius  was  fl/e2)  was  0.8  mm. 


Spatial  and  Thermal  Resolution  of  2D  MRT 

Measuring  two-dimensional  temperature  gradients  in  the  layered  phantom  during  infrared  laser 
exposure  assessed  the  spatial  resolution  of  the  MRT  method.  A  difference  in  the  spatial 
temperature  gradients  due  to  the  non-homogeneity  of  the  melanin  layer  in  the  phantom  was  readily 
apparent.  Infrared  imaging  revealed  the  optical  propagation  of  the  laser  beam  through  the  phantom. 
The  beam  propagated  with  a  small  amount  of  scattering  and  very  little  absorption  through  the  clear 
agarose  gel,  but  was  highly  scattered  and  absorbed  in  the  melanin  layer  (Figure  3A).  A 
corresponding  MR  magnitude  image  of  the  phantom  is  shown  for  comparison  in  Figure  3B. 


47 


Ct**r*9arot*g«( 

Layer  of  agaro**  and 
RPE  rmUnotomn 

«W  nm  t***r  btam 


Figure  3.  (A)  Infrared  image  taken  of  the  806  nm  laser  beam  propagating  through  the  layered 
phantom.  The  IR  image  was  taken  with  a  Sony  F717  camera  (IR  imaging  mode)  of  the  phantom 
outside  of  the  MR  cavity.  During  the  actual  experiment,  the  laser  beam  was  normally  incident  on  the 
melanin  layer.  (B)  MR  magnitude  image  of  the  layered  phantom.  The  melanin  layer  appears  brighter 
than  the  surrounding  clear  agarose  because  of  melanin’s  inherent  T1  contrast. 


Because  of  its  paramagnetic  moment,23  the  melanin  layer  has  higher  T1  contrast  than  the 
surrounding  agarose.  Temperature  gradients  were  measured  in  the  layered  phantom  using  MRT 
during  laser  exposures  delivering  sufficient  power  into  the  phantom  to  produce  heating  (Figure  4). 
A  thermographic  image  was  constructed  by  obtaining  the  phase  difference  between  a  magnetic 
resonance  image  taken  prior  to  the  onset  of  the  laser,  and  one  taken  during  the  laser  exposure.  The 
thermographic  (phase  difference)  image  is  shown  in  Figure  4A.  The  spatial  resolution  of  the  MRT 
images  was  determined  by  the  number  of  image  voxels  within  the  physical  field  of  view.  The  in¬ 
plane  resolution  of  the  2D  MR  images  was  0.25  mm.  The  MRT  spatial  resolution  was  sufficient  to 
distinguish  between  temperature  gradients  measured  within  the  plane  of  the  melanin  layer,  and 
those  orthogonal  to  the  layer,  as  shown  in  Figure  4B.  The  gradient  within  the  melanin  layer  is  wider 
compared  to  the  orthogonal  profile.  That  difference  can  be  explained  by  the  low  optical  absorption 
of  the  laser  beam  within  the  clear  agarose,  in  contrast  to  the  high  optical  absorption  -  and  hence, 
heating  -  of  the  laser  beam  in  the  melanin  layer. 

The  theoretical  temperature  resolution  of  MRT  based  on  the  PRF  method  may  approach  0.1  °C.12 
By  using  an  external  water  bath  to  heat  the  phantom  uniformly  to  a  known  temperature,  the  MRT 
measurements  were  calibrated  so  that  an  absolute  temperature  could  be  derived  from  the 
temperature-dependent  phase  changes  in  the  image.  In  addition,  by  comparing  the  agreement  of 
MRT  of  the  water-bath  heated  phantom,  determined  from  2D  images,  with  temperature 
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measurements  made  with  a  Vitek  probe,  an  estimate  of  the  accuracy  of  the  method  was  attained, 
and  was  found  to  be  within  0.5  °C,  at  least  for  temperature  excursions  less  than  10  °C. 


Figure  4.  Orthogonal  temperature  profiles  in  the  layered  phantom  during  laser  exposure.  (A)  The 
phase  image  of  the  phantom  during  laer-induced  heating:  the  dark  band  near  the  top  of  the  image  is  an 
artifact  in  the  phase  unwrapping  algorithm.  The  white  lines  in  the  figure  indicate  the  relative 
orientations  of  the  temperature  gradients  with  respect  to  the  melanin  layer  in  the  phantom.  (B)  A 
graph  of  the  temperature  gradients  in  the  phantom  as  a  function  of  displacement  from  the  focal  point 
of  the  laser  beam  in  the  melanin  layer.  The  filled  symbols  represent  MRT  measurement  in  a  slice 
extending  through  the  hottest  point  in  the  phantom,  i.e.  within  the  plane  of  the  melanin  layer  through 
the  laser  focal  point;  open  symbols  represent  MRT  measurements  in  a  slice  running  through  the  laser 
fecal  point,  in  c.  plane-perpendicular  to  the  melanin  Inysr.  Both  dices  extended  across  the  entire 
width  of  the  phantom.  Continuous  lines  through  both  data  sets  are  best-tit  Gaussian  functions.  The 
temperature  gradient  within  the  melanin  layer  is  broader  than  the  gradient  across  the  layer  extending 
into  the  clear  agarose  on  each  side  of  the  melanin  layer,  due  to  the  greater  scattering  and  absorption  of 
the  laser  energy  in  the  melanin  granules. 


Three-Dimensional  Temperature  Profiles 


Once  at  the  thermal  equilibrium,  volume  image  scans  were  made  to  construct  a  three-dimensional 
thermographic  image  of  steady-state  heat  flow  through  the  phantom,  starting  from  the  focal  point  of 
the  laser  beam  where  it  interacted  with  the  melanin  layer.  The  spatial  resolution  of  the  3D  images 


was  0.25  mm  in  all  three  axes.  The  data  sets  shown  in  Figure  5  illustrate  the  temperature 
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topography  inside  the  volume  of  the  phantom.  It  may  be  appreciated  that  the  temperature  gradients 
follow  a  Gaussian  profile  due  to  the  Gaussian  distribution  of  power  in  the  laser  beam.  The  diameter 
of  the  laser  beam  irradiating  the  phantom  also  affected  the  shape  of  the  three-dimensional  heating 
profile.  The  three-dimensional  thermogram  resulting  from  exposure  to  the  0.8  mm  radius  laser 
beam,  shown  in  Figure  5A,  is  narrower  than  the  thermogram  produced  by  the  2.4  mm  radius  beam, 
shown  in  Figure  5B.  The  SNR  for  the  three-dimensional  magnitude  image  was  11.5  (calculated  as 
the  inverse  of  the  standard  deviation  of  the  phase  over  the  ROI).  To  obtain  an  estimate  of  the  error 
in  the  temperature  measurements,  the  squares  of  the  errors  across  the  image  ROI  were  added,  and 
the  square  root  taken.  Converting  back  to  temperature,  the  derived  estimate  of  error  in  the  3D  data 
set  was  ±1 .9  degrees. 


Figure  5.  Surface  plots  taken  from  three-dimensional  maps  of  the  layered  phantom  during  steady-state 
laser  exposure.  (A)  The  thermogram  resulting  from  exposure  of  the  phantom  to  the  0.8  mm  radius  laser 
beam.  (B)  The  thermogram  resulting  from  exposure  of  the  phantom  to  the  2.4  mm  radius  laser  beam. 
The  inset  in  the  upper  right-hand  comer  of  each  figure  shows  the  phase  difference  image  (thermogram) 
of  the  melanin  layer  (the  image  slice  was  taken  within  the  plane  of  the  melanin  layer).  The 
thermograms  illustrate  the  temperature  gradients  spreading  out  into  the  phantom,  away  from  the  center 
of  the  laser  beam  incident  on  the  melanin  layer.  The  view  of  the  surface  plot  is  approximately  90°  to 
the  axis  of  the  incident  laser  beam. 


Correlation  of  Thermal  Model  and  MRT  Measurements 

In  order  to  compare  a  temperature  profile  measured  by  MRT  during  laser  exposure  with  the  heating 
gradient  predicted  by  the  heat  diffusion  model  described  in  Section  2.4,  a  two-dimensional  image 
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slice  was  taken  out  of  the  three-dimensional  data  set  shown  in  Figure  5A.  The  temperature  profile 
obtained  from  this  image  slice  is  shown  in  Figure  6  as  the  discrete  points,  while  the  continuous  line 
in  Figure  6  shows  the  temperatures  predicted  by  equation  (5).  The  error  bars  shown  in  Figure  6 
represent  the  error  estimates  derived  from  the  three-dimensional  data  sets  as  described  above 
(“Three-Dimensional  Temperature  Profiles”).  The  phantom  heating  was  produced  by  continuous 
wave  laser  exposure  of  21.6  W/cm2,  0.8  mm  beam  radius  (1/e2),  and  700  seconds  duration.  The 
maximum  temperature  rise  at  the  center  of  the  phantom  was  about  32  °C,  which  was  extreme  for 
hyperthermia  treatment,  but  a  good  test  of  the  ability  of  MRT  to  record  a  large  temperature 
excursion  with  high  resolution.  The  agreement  in  Figure  6  between  the  theoretical  and 
experimental  values  is  quite  close.  This  result  indicates  that  the  observed  values  were  well 
correlated  with  the  predicted  values,  and  that  the  heat  diffusion  model  is  an  excellent  predictor  of 
heat  transfer  in  such  basic  applications. 


Figure  6.  Comparison  of  the  heat  diffusion  model  computational  results  and  the  temperature  profile 
measured  by  MRT  within  the  layered  ocular  phantom  during  IR  laser  exposure.  The  filled,  discrete 
symbols  represent  the  MRT  measurements  taken  from  a  central  image  slice  of  the  3D  data  set  shown 
in  Figure  5A.  The  continuous  line  represents  the  temperature  rise  in  the  phantom  predicted  by  the 
thermal  model  as  a  function  of  displacement  from  the  center  of  the  laser  beam  incident  on  the  melanin 
layer.  Error  bars  on  the  discrete  points  are  estimates  of  temperature  measurement  standard  deviation, 
calculated  as  described  in  the  section,  titled  “MR  Thermography”. 
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DISCUSSION 


Practical  Spatial,  Temperature,  and  Temporal  Resolution  of  MRT 

The  accuracy  of  the  temperature  calculations  made  with  MRT  depend  on  several  factors,  the  most 
important  of  which  is  the  determination  of  the  thermal  coefficient,  a.  In  the  present  application,  for 
elevations  of  less  than  10  °C  produced  by  passive  heating  of  the  uniform  phantom  with  the  water 
bath,  the  temperature  values  determined  by  MRT  were  compared  to  measurements  made  with  a 
Vitek  fiber  optic  probe  having  an  accuracy  of  0.1  °C,  and  usually  agreed  to  within  0.5  °C.  The 
temperature  resolution,  however,  may  depend  on  the  magnitude  of  the  temperature  excursion.  For 
example,  during  the  more  extreme  heating  produced  by  the  laser  exposure,  the  SNR  of  the 
temperature  measurements  decreased,  and  for  temperature  changes  of  more  than  20  °C,  the 
resolution  decreased  in  2D  MRT  to  about  1.4  °C,  and  to  1.9  °C  for  3D  MRT.  This  implies  that  the 
practical  temperature  resolution  of  the  MRT  method  decreases  with  an  increase  in  the  relative 
temperature  change  in  the  target,  although  the  resolution  is  probably  better  than  1  °C  for  clinically 
relevant  temperature  excursions.  Another  potential  source  of  error  is  due  to  an  interaction  between 
the  size  and  orientation  of  the  heat  source  with  the  external  magnetic  field.  Peters  et  al. 
demonstrated  that,  with  an  interstitial  laser  fiber  having  a  1-cm  long  quartz  diffuser,  there  was  an 
orientation-dependent  effect  on  the  measured  temperature  that  was  ascribed  to  a  change  in  the 
volume-magnetic  susceptibility  of  the  tissue  along  the  axis  of  the  fiber  diffuser. 13,24,25  With  a  laser 
beam  applied  externally,  as  in  the  present  case,  this  orientation  artifact  would  be  expected  to  be 
minor,  because  the  effective  heat  source  in  the  phantom,  i.e.  the  image  of  the  laser  beam,  is  thin 
compared  to  the  optical  zone  absorbing  the  laser  energy,  i.e.  the  melanin  layer.  In  a  more 
complicated  target  such  as  the  eye  (in  which  the  RPE  layer  containing  melanin  is  very  thin),  a 
correction  for  orientation  dependence  may  have  to  be  introduced.  This  will  be  examined  in  future 
work. 

In  terms  of  spatial  resolution,  an  effective  resolution  of  0.25  mm  was  obtained.  Given  the  average 
thickness  of  the  human  retina  (350-500  microns),  the  RPE  (20  microns),  and  the  choroid  (100-200 
microns),26  the  0.25-mm  spatial  resolution  in  the  present  MRT  implementation  is  sufficient  to 
provide  two  or  three  temperature  measurement  points  through  the  combined  thickness  of  the  retina 
and  choroid  layers,  although  insufficient  to  distinguish  retinal  and  RPE  temperature.  Nevertheless, 
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the  currently  achievable  resolution  supports  the  application  of  MRT  to  monitor  the  heating  of  the 
target  vascular  membrane,  as  well  as  surrounding  retina  and  RPE,  during  treatment  of  choroidal 
neovascularization  by  TTT.1'3  A  limitation  of  TTT  at  present  is  that  the  target  temperature  required 
for  the  optimal  therapeutic  effect  has  not  been  established.  Use  of  MRT  for  feedback  control  during 
TTT  could  provide  critical  temperature  information  to  correlate  with  the  resulting  tissue  responses 
to  the  treatment. 

Validation  of  Thermal  Models 

Experiments  in  this  investigation  indicated  a  good  agreement  between  the  calculated  data  from  the 
basic  thermal  model  and  data  measured  with  the  MRI  technique.  Such  an  agreement  depends 
critically  on  the  physical  and  optical  properties,  such  as  the  thermal  conductivity,  specific  heat 
capacity,  density,  absorption  and  scattering  coefficients,  and  the  pigmentation  and  the  homogeneity 
of  the  tissue,  that  are  entered  into  the  model  and  control  the  calculation  of  both  the  radiation 
distribution  and  the  thermal  energy  distribution.  The  near  infrared  laser  has  been  reported  to  have 
less  scattering  effect  in  pigmented  tissue  compared  to  the  visible  wavelengths,  although  the 
scattering  is  still  the  dominant  process  at  806  nm.  The  effective  penetration  depth  at  806  nm  in  the 
melanized  layer  of  our  phantom  has  been  calculated  from  equation  (5)  to  be  ~0.7  mm,  which  is  less 
than  the  width  of  the  phantom’s  melanin  layer  (2  mm  thick).  Consequently  most  of  the  optical 
heating  by  the  laser  radiation  at  806  nm  was  confined  in  the  slab  of  the  melanin  layer,  consistent 
with  the  light  propagation  revealed  with  the  infrared  image  shown  in  Figure  3 A.  This  also  will  slow 
the  laser-induced  heating  time  in  the  thick  melanin  layer  of  the  phantom  compared  to  the  thin  RPE 
layer  (~10  pm)  of  the  physiological  eye.  Another  factor  slowing  the  heating  process  in  the  melanin 
layer  is  the  size  of  the  phantom.  Within  the  2.5cm  radius  of  the  phantom,  more  heat  will  be 
distributed  by  diffusion  throughout  the  agarose  volume  surrounding  the  melanin  layer.  But  that  is 
analogous  to  the  real  eye,  where  heating  the  ocular  tissue  also  involves  heating  the  vitreous  body. 

Variation  in  the  concentration  of  melanin  granules  in  the  absorbing  layer  will  cause  variation  in  the 
optical  properties,  which  can  critically  affect  the  thermal  model’s  ability  to  predict  temperature 
distribution.  To  date,  there  have  been  few  consistent  measurements  of  optical  properties  of  the 
retinal  pigment  epithelium  (RPE)  or  choroid  layer.27  Accurate  measurements  of  the  optical 
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properties  are  critical  inputs  for  heating  models,  and  better  validation  of  thermal  models  can  be 
done  by  comparing  theoretical  predictions  of  infrared  laser  heating  with  MRT  measurements  using 
different  concentration  of  melanin  granules  in  absorbing  layers  of  phantom  models.  Another 
consideration  is  that  the  influence  of  radiational  transport  will  be  less  important  to  the  treated 
volume  in  the  case  of  long  heating  time  (because  of  thermal  equilibration). 

The  present  model  is  conceptually  straight  forward,  with  further  computational  frugality  achieved 
by  taking  simple  considerations  of  boundary  conditions.  The  geometric  and  thermal  assumptions 
reflect  values  which  limit  the  accuracy  and  flexibility  of  the  model.  More  complex  models  can  be 
developed  for  future  laser  hyperthermia  studies  using  a  finite  element  approach.  An  advanced 
model  would  take  into  consideration  radiational  transport  in  ocular  tissue  among  different  layers, 
such  as  the  RPE  and  choroid,  with  different  optical  and  physical  properties.  In  addition  blood  flow 
in  the  choroid  would  be  another  factor  to  consider  in  in  vivo  studies. 

Technical  Limitations  and  Improvements  in  the  Method 

The  temperature  and  spatial  resolution  of  the  present  system  appear  to  be  adequate  for  application  to 
the  human  eye.  The  temporal  resolution,  however,  could  be  improved.  There  is  a  tradeoff  between 
temporal,  spatial,  and  temperature  resolution.  Increasing  the  resolution  of  one  inevitably  reduces 
the  resolution  of  one  or  both  of  the  remaining  parameters;  therefore,  the  experiment  must  be 
designed  to  maximize  the  most  important  parameter.  For  example,  in  a  study  of  experimental 
interstitial  laser  hyperthermia  therapy  to  the  prostate,  a  single  image  slice  was  acquired  in  6.7  sec, 
but  with  a  slice  thickness  of  5  mm.28  Certainiy,  a  faster  image  acquisition  rate  would  be  very  useful 
in  improving  the  accuracy  of  heating  time  course  measurements  in  the  target  tissue.  It  is  likely  that 
advances  in  MR  instrumentation  will  improve  overall  imaging  performance. 

Application  of  MRT  to  Characterizing  the  Interaction  of  IR  Laser  with  Ocular  Tissue  Components 
The  MR  measurements  in  the  ocular  phantom  provided  a  means  for  directly  studying  the  factors 
governing  the  temperature  distribution  during  infrared  laser  exposure.  The  physical  distribution  of 
the  melanin  pigment  granules,  the  concentration  (or  density)  of  melanosomes  in  the  tissue,  and  the 
thickness  of  the  pigmented  tissue  layer  all  determined  the  effective  optical  zone  in  which  laser- 
tissue  interaction  occurred.  There  were  limitations  in  the  ability  of  the  phantom  to  simulate  thermal 
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interactions  likely  to  occur  in  the  living  eye.  Practical  considerations  in  the  construction  of  the 
layered  phantom  limited  its  ability  to  simulate  the  eye  morphology  more  accurately.  For  example, 
the  laser-induced  heating  time  course  measured  in  these  experiments  was  considerably  slower  than 
the  thermal  behavior  predicted  for  retina  and  RPE  during  laser  irradiation.29,30  Considering  that  the 
heat  diffusion  process  will  be  complicated  for  a  blood-perfused,  inhomogeneous  organ  such  as  the 
human  eye,  MRT  provides  a  non-invasive  means  for  determining  in  vivo  the  heating  time  course  as 
well  as  the  equilibrium  temperature,  both  essential  for  understanding  the  biological  response  of  the 
tissue  to  thermal  stress. 

The  interaction  of  the  806  nm  CW  diode  laser  with  melanin  is  dominated  by  thermal  effects,  and 
depends  primarily  on  photon  scattering.  For  example,  the  relatively  high  absorption  coefficient  of 
melanin  produced  a  high  temperature  rise  in  the  small  volume  of  the  melanin  layer,  which  caused 
heat  to  diffuse  throughout  the  phantom.  This  effect  is  illustrated  in  Figure  4,  which  shows  that  the 
laser  radiation  was  mostly  absorbed  within  the  melanin  layer,  and  that  heating  profiles  through  the 
phantom  depended  on  the  particular  geometry  of  the  absorbing  layer  in  the  phantom.  Selection  of 
suitable  laser  parameters  for  hyperthermia  therapy,  therefore,  will  depend  on  the  ability  to  predict 
accurately  the  interaction  of  the  laser  beam  with  the  tissue  and  the  resultant  heat  flows  through  the 
treated  and  adjacent  non-treated  tissue.  Such  predictions  will  require  the  use  of  a  validated 
treatment  model.  The  good  agreement  of  the  MRT  measurements  with  the  predictions  of  a  classical 
heat  diffusion  model  show  that  an  important  use  of  the  method  is  to  validate  theoretical  models  of 
laser -tissue  thermal  interaction.  One  practical  embodiment  of  this  ability  would  be  to  support  the 
clinical  use  of  a  thermal  model  to  plan  therapy  and  optimize  the  exposure  parameters  for  laser- 
induced  hyperthermia  treatment  protocols. 
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